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ABSTRACT 


An  evaluation  of  the  currently  accepted  visco- 
inertial  flow  theory  has  indicated  that  in  normal  laboratory 
determination  of  the  absolute  permeability  and  the  inertial 
resistance  coefficient  a  correction  for  molecular  streaming 
is  warranted.  Modified  forms  of  the  visco-inertial  theory 
have  been  developed  considering  molecular  streaming  for  flow 
under  both  low  and  high  pressure  gradients. 

A  graphical  means  of  solution  for  the  aforementioned 
parameters  of  a  medium  has  been  proposed  for  flow  at  a  low 
pressure  gradient,  and  a  general  numerical  technique  has  been 
developed  for  use  with  both  modified  flow  equations. 

It  has  been  shown  that  the  theory  developed  also 
applies  for  gas  flowing  in  the  presence  of  a  static  liquid 
phase.  Experimental  determination  of  the  effective  perme¬ 
ability,  the  Klinkenberg  b-factor  and  the  coefficient  of 
inertial  resistance  has  indicated  that  these  parameters  are 
continuous  functions  of  the  liquid  saturation.  While  the 
effective  permeability  has  been  found  to  decrease  with  satura¬ 
tion  and  the  Klinkenberg  b-factor  to  remain  relatively  con¬ 
stant,  the  inertial  resistance  coefficient  has  been  found  to 
increase  at  an  increasing  rate. 
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INTRODUCTION 


For  some  time  now  flow  of  a  single  phase  fluid 

through  saturated  porous  media  has  been  studied  rigorously. 

It  has  been  customary  to  describe  flow  in  the  viscous  and 

visco-inertial  flow  regions  by  means  of  Darcy's  Law  and  the 

(1  2 ) 

Forchheimer  equation  respectively  '  .  The  Klinkenberg 

b-f actor,  the  permeability  and  the  inertial  resistance  co¬ 
efficient  are  empirical  constants  used  to  characterize  a 
porous  medium. 

Generally  the  presence  of  a  connate  water  is  in¬ 
herent  to  those  formations  which  contain  petroleum  hydro¬ 
carbon  deposits.  Hence,  it  is  important  to  predict  the  flow 
behavior  of  a  gas  in  the  presence  of  a  liquid  such  as  a  con¬ 
nate  water  or  a  hydrocarbon. 

The  theory  describing  flow  in  a  dry  porous  medium 
has  been  extended  to  describe  the  phenomenon  of  flow  in  the 

presence  of  a  second  phase.  Under  viscous  flow  conditions 

( 3 ) 

this  type  of  analysis  has  been  verified  experimentally 

The  conclusions  which  have  been  made  are  that  the  permeability 

to  the  gas  is  a  continuously  decreasing  function  of  the  liquid 

saturation,  unique  for  each  sample,  while  the  Klinkenberg 

b-factor  remains  essentially  constant.  An  attempt  has  also 

been  made  to  apply  the  visco-inertial  flow  theory  to  gas 

(4) 

flow  in  the  presence  of  a  second  phase 


2 


The  object  of  this  study  was  to  theoretically  and 
experimentally  verify  the  visco-inertial  flow  theory  and  to 
determine  the  effect  that  an  increasing  liquid  saturation 
has  upon  the  characteristic  parameters  of  a  porous  medium 
under  visco-inertial  flow  conditions. 


3 


LITERATURE  REVIEW 


Visco-Inertial  Flow 

( 5 ) 

According  to  Scheidegger  ,  Forchheimer,  drawing 
an  analogy  with  the  phenomenon  occurring  in  tubes,  sug¬ 
gested  that  high  velocity  flow  through  porous  media  may  be 
described  by  an  expression  of  the  form 


aq  +  bq 


2 


(1) 


where  q  is  the  seepage  velocity  and  a  and  b  are  thought  to 
be  constants. 

From  a  development  based  on  Reynold's  theory  of 
resistance  to  flow  and  the  Kozeny  assumption  that  a  porous 
medium  is  analogous  to  parallel  capillaries  with  common 
hydraulic  radius,  Ergun  and  Orning v  }  presented  an  expression 
of  the  form 


cyq  +  dpq 


2 


(2) 


where  c  and  d  are  constants  and  functions  of  the  medium. 
They  illustrated  the  validity  of  this  expression  by  showing 
the  linearity  of  (p^  -  versus  pqQ  for  flow  through 

unconsolidated  porous  media  of  fluids  which  were  assumed  to 
behave  ideally. 

( 7 ) 

Green  and  Duwez  characterized  the  structure  of 
a  porous  medium  with  two  length  parameters  F^  and  Fg ,  using 


4 


the  quadratic  equation 

dp 

"  —  =  FAyq  +  FBpq  (3) 

dx 

is  the  viscous  resistance  coefficient  and  Fg  is  the  iner¬ 
tial  resistance  coefficient ,  both  independent  of  the  fluid 
properties.  Using  a  form  of  Equation  (3)  integrated  for  a 
constant  mass  flow  rate  of  an  ideal  gas,  they  evaluated  the 
resistance  coefficients  for  four  sintered  porous  metal  speci¬ 
mens  . 

( 8 ) 

Cornell  and  Katz  v  conducted  a  study  of  "inertial, 
quasi  turbulent  or  turbulent"  flow  of  various  gases  through 
twenty-four  consolidated  samples  of  sandstone,  limestone, 
and  dolomite.  They  concluded  that  an  expression  of  the  form 
of  Equation  (3)  is  valid. 

The  motion  of  ordinary  Newtonian  fluids  can  be 

described  by  means  of  the  Navier-Stokes  equation.  In  a 

porous  medium  gross  simplifications  must  be  made  in  order  to 

( 9 ) 

formulate  boundary  conditions.  Irmay  analyzed  an  idealized 
pore  system  of  spherical  particles,  and  Tek^^  represented 
the  flow  channels  by  a  sequence  of  truncated  cones,  giving 
rise  to  successive  restrictive  orifices  along  the  channel. 

They  both  indicated  the  validity  of  an  expression  similar  to 
Equation  (3)  which  may  be  written  as 

-  ^  =  Famc  +  FBpq|q|  (4) 


. 


5 


where  F,  and  F_  are  dependent  upon  the  nature  of  the  pore 
geometry.  The  conditions  imposed  in  arriving  an  Equation  '4) 
are  horizontal  steady  linear  isothermal  flow  of  a  single 
phase  Newtonian  fluid  through  chemically  inactive  homogeneous 
isotropic  geometrically  stable  saturated  media. 

Equation  (3) ,  the  differential  form  which  will  be 
referred  to  as  the  Forchheimer  equation,  has  been  put  into  a 
more  useful  form,  by  combining  with  the  modified  gas  law  and 
integrating  for  a  constant  mass  flow  rate.  Realizing  that 
viscosity  y  and  compressibility  Z  are  pressure  dependent  under 


isothermal  conditions,  Mackett 


(11) 


used 


Z  T  d  A 
O  O  "  c 


D 

-  r-.  p 

(  /  1  ^ 


p  y  T  Q  T  L  o 
*0  o  c  o  r 


Zy 


dD  - 
-  r 


P-r 

'2  Pr 


dD  ) 

*■  -v 


F  +  F 
~  A  B 


Q  p 
o  o 


Ay  . 
ia 


(5) 


where  the  integrated  average  viscosity  was  taken  as 


Pr 


."1  Pr  fT2  Pr 

—  dp  -  —  dp 


*  r 


y .  =  y  - 

la  o  pr 


Pr. 


(6) 


;Ai  pr  %  2  pr 

J  Zy  ^r  J  Zy  - r 
o  r  o  r 


y 

where  y  =  — 

r 

1 1 

~  o 

He  evaluated  the  isothermal  variation  of  viscosity 
and  compressibility  of  nitrogen  with  pressure  to  investigate 
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the  magnitude  of  errors  resulting  from  the  assumption  of 

average  fluid  properties.  He  concluded  that  under  normal 

laboratory  conditions  it  was  valid  to  use  values  of  viscosity 

and  compressibility  evaluated  at  the  arithmetic  mean  pressure 

and  temperature.  This  indicated  that  the  integrated  form 

( 12 ) 

used  by  Cornell  is  valid 

M(p,  2  -  P92)  w 

_  -  =  F  +  F  —  (7) 

2y  ZRTL  w/A  Ay 


The  non-rigorous  nature  of  the  derivation  of  the 
Forchheimer  equation  left  open  the  question  whether  the 
factors  F  and  FD  were  geometrical  and  structural  variables 
or  were  dependent  upon  pressure  or  velocity.  Since  structure 
itself  may  change  to  a  certain  extent  with  varying  ambient 
conditions,  changes  in  temperature  and  pressure  may  affect 


fa  and  FB- 


(13) 

Greenberg  and  Weger  studied  the  constancy  of 

F  and  F  for  flow  through  porous  metals  at  high  rates  and 
pressures  up  to  2000  psi.  They  concluded  that  these  coeffi¬ 
cients  remained  constant  under  isothermal  conditions. 


Viscous  Flow 

The  criterion  for  kinematic  similarity  has  been 

( i4 ) 

applied  by  Hubbert  to  define  the  Reynolds  Number,  NRe , 

which  is  a  measure  of  the  ratio  of  inertial  effects  to  the 
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viscous  effects. 


From  Equation  (7)  this  may  be  written  as 


N 


Re 


fa  ^ 


(8) 


At  very  low  rates,  for  example  at  Reynolds  Numbers 
less  than  0.01,  the  portion  of  the  pressure  drop  due  to 
inertial  effects  is  very  small  compared  to  that  due  to  the 
viscous  effects.  Hence,  the  inertial  effects  may  be  neglected 
in  Equation  (3)  resulting  in 

dp 

- -  F  yq  (9) 

dx 

(15) 

When  a  comparison  is  made  with  Darcy's  Law 

dp  1 

-  —  =  -  yq  (10) 

dx  k 


the  viscous  resistance  coefficient  is  found  to  be  the  inverse 
of  permeability  k. 

For  linear  flow  of  a  constant  mass  flow  rate, 

(16) 

Mackett v  '  again  indicated  the  validity  of  using  values  of 
the  viscosity  and  compressibility  evaluated  at  the  mean  tem¬ 
perature  and  pressure.  The  resulting  integrated  forms  are 

m(p-l2  -  p22) 

2yZRTL  w/A 

AToZo(Pl2  -  P22) 

2^7z  t  po  lqo 


1 

k 


(11) 


k 


(12) 
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Following  the  investigation  of  Kundt  and  Warburg, 

( 17 ) 

Klmkenberg  observed  that  as  the  mean  free  path  length 

of  the  fluid  approaches  the  dimensions  of  the  flow  conduit, 
the  gas  layer  immediately  adjacent  to  the  wall  of  a  capillary 
has  a  finite  velocity  in  the  direction  of  flow.  The  mean  free 
path  length  or  the  average  distance  travelled  by  a  gas  mole¬ 
cule  between  intermolecular  collisions  is  a  function  of  the 

(18) 

pressure,  temperature  and  the  nature  of  the  gas v  . 

1  RT 

X  =  - 5-  =  -3 - T  (13) 

/2  ird  n  / 2  upNd 


Drawing  an  analogy  to  the  flow  in  glass  capillary 

tubes,  Klinkenberg  corrected  the  apparent  permeability  of  a 

porous  medium  to  the  absolute  permeability  exhibited  to  a 

fluid  of  a  very  low  mean  free  path.  On  the  basis  of  the 

theory  and  extensive  experimentation  he  proposed  that  the 

apparent  permeability  k  ,  calculated  from  Darcy's  law  using 

a 

Equations  (11)  and  (12) ,  be  related  to  the  absolute  perme¬ 
ability  by 

4cA 

k  =  k  (1  +  - )  (14) 

cl 

r 

For  any  gas  c  is  a  constant  and  is  approximately 
equal  to  1.0.  Substituting  for  the  mean  free  path  length 
in  Equation  (14) 


4cRT 

/2  TTd2Npr 


k  (1  + 


(15) 
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4  CRT 

Under  isothermal  conditions  the  factor  — — - x — 

/2  tt d  Nr 

is  a  constant  for  a  given  gas  and  medium.  This  is  referred 
to  as  the  Klinkenberg  b-f actor,  which  is  inversely  proportion¬ 
al  to  an  effective  path  radius  r.  Therefore  equation  (15) 
may  be  written 

b 

k  =  k(l  +  -)  (16) 

a.  — 

p 

where  p  is  the  arithmetic  mean  pressure. 

Emphasis  has  been  placed  on  slippage  because  of 

its  presence  under  normal  laboratory  conditions.  Stewart  and 
(19) 

Owens  performed  tests  on  heterogeneous  limestones  to  in¬ 

vestigate  the  extent  to  which  slippage  and  inertial  effects 
can  exist  simultaneously.  Slippage  was  controlled  by  main¬ 
taining  constant  mean  free  path  conditions.  At  very  large 
Reynolds  Numbers  the  slippage  effect  disappeared.  Neverthe¬ 
less  at  high  flow  rates  it  was  found  that  both  slippage  and 
inertial  effects  contributed  significantly  to  the  phenomenon 
observed . 

Gas  Flow  in  the  Presence  of  a  Second  Phase 

Experimental  evidence  has  indicated  that  relationships 
similar  to  those  obtained  to  predict  gas  flow  behavior  through 
dry  porous  media  are  valid  in  describing  the  flow  of  gas  in 
the  presence  of  a  liquid  phase. 
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Rose^^  and  Estes  and  Fulton  showed  that 

Klinkenberg ' s  interpretation  of  slippage  could  be  applied 
to  both  homogeneous  gas  flow  and  flow  of  gas  in  a  liquid 
saturated  porous  medium.  It  was  observed  that  under  iso¬ 
thermal  conditions  the  relative  effect  of  slippage  was  con¬ 
stant  as  the  Klinkenberg  b-factor  remained  relatively  con¬ 
stant  for  various  saturations.  The  permeability  was  found 

to  degrease  continuously  with  increasing  liquid  saturation. 

(22) 

Hamm  and  Eilerts  studied  the  effects  of  a  con¬ 

densate  saturation  in  typical  sandstone  and  limestone  reser¬ 
voir  rocks  on  the  mobility  of  gas.  Under  viscous  flow  con¬ 
ditions  they  investigated  the  effects  of  pressure,  liquid 
saturation  and  velocity  on  the  mobility,  which  is  defined  as 
the  ratio  of  the  permeability  to  the  viscosity.  They  con¬ 
cluded  that  increased  gas  velocity,  increased  liquid  satura¬ 
tion  and  increased  pressure  would  decrease  the  mobility. 
These  effects  are  attributable  respectively  to  inertial  ef¬ 
fects,  the  influence  of  the  liquid  saturation  and  the  reduc¬ 
tion  of  slippage  effects  and  increase  of  gas  viscosity  with 
increasing  pressure. 

Using  a  sample  of  heterogeneous  limestone  which 

(23) 

exhibited  negligible  slippage  effects,  Stewart  and  Owens 
studied  two-phase  flow  of  gas  and  oil  under  external  gas 
drive.  They  found  two  regions  of  flow  similar  to  those 
which  are  evident  in  single-phase  flow. 


r>  :ji  a  nox3  r.i  -? 

. 

. 

-siujfia  bxu  xl  b  r  .  uni  ,y^xdoI9V  asp  baasD^oni  ^s  ^  b  >u.  o 
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Iffly  and  Naville 


(24) 


presented  an  equation,  describ¬ 


ing  visco-inertial  flow,  which  was  accredited  to  Houpeurt 


v  ( 1  +  <f>U  —  I  v  j  )  = 

V 


-  —  grad  p 

<P  y 


(17) 


U  and  k^,  the  effective  permeability  to  the  gas,  are  the  two 
length  parameters.  The  integrated  form  of  the  equation  for 
linear  flow  of  an  ideal  gas  was  used 


U  w  k  A  ^2  k  A  M(p,^-p  2) 

w(l  +  - )  =  -a-  1  pdp  =  -2 - li-L* -  (18) 

y  A  y  p^  2y  LRT 


Flow  tests  were  performed  on  a  sandstone  sample  at  various 

saturations  of  a  frozen  immobilized  second  phase.  For  all 

2  2 

saturations  a  plot  of  (p^  -p2  ) /w  versus  w  appeared  to  be 
linear.  The  parameters  U  and  kg  were  evaluated  for  all 
saturations.  With  an  increasing  saturation  of  the  second 
phase  kg  was  found  to  decrease  continuously  while  U  increased 
in  magnitude. 

Results  were  found  to  compare  favorably  with  those 

obtained  for  simultaneous  flow  of  two  phases.  This  procedure 

was  not  recommended  because  of  instability,  but  served  to 

indicate  that  the  phenomenon  observed  occurs  independently  of 

the  saturating  medium. 

(25) 

Nichol  also  studied  visco-inertial  flow  behavior 

in  the  presence  of  an  immobile  second  phase.  He  investigated 


12 


the  variation  of  the  inertial  resistance  coefficient  F  ,  as 

B 

defined  in  Equation  (7) ,  with  an  increasing  saturation  of 
glycerin,  the  second  phase.  The  coefficient  was  found  to 
decrease  slightly  at  low  saturations  but  increase  sharply 
thereafter . 
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EVALUATION  OF  THE  CURRENT  THEORY 

Analysis  of  Laboratory  Flow  Tests 

The  dynamic  flow  equations  may  be  applied  to  predict 
gas  flow  behavior  in  a  porous  medium  under  the  conditions  of 
viscous  and  visco-inertial  flow  provided  the  three  parameters 
absolute  permeability,  the  Klinkenberg  b-f actor  and  the  in¬ 
ertial  resistance  coefficient  are  known.  It  has  been  found 
advantageous  to  obtain  these  parameters  from  laboratory  tests 
upon  representative  samples. 

At  first  the  permeability  of  a  test  specimen  was 
determined  by  conducting  a  single  laboratory  flow  test  and 

applying  Darcy's  Law  to  the  resulting  data. 

(26) 

Klinkenberg k  ;  observed  that  this  apparent  per¬ 
meability  of  a  porous  medium  to  a  gas  was  dependent  upon 
the  nature  of  the  medium,  the  mean  flowing  pressure  and  the 
nature  of  the  gas.  He  proposed  that  the  absolute  permeability 
be  determined  by  means  of  a  multiple  point  flow  test.  By 

Equation  (16) ,  from  a  plot  of  k  versus  1/p,  termed  a  Klin- 

a 

kenberg  plot,  the  absolute  permeability  is  the  extrapolated 
intercept  and  the  Klinkenberg  b-factor  is  the  slope  divided 
by  the  absolute  permeability. 

For  flow  in  the  visco-inertial  region  Equation  (7) 
may  be  modified,  recalling  that  is  the  reciprocal  per¬ 
meability. 


I3T  jo  aHT  go  mc  : f:  i  A'_- 


;o  a  j  j±sc.  . 


I  3  anoxdxbnoo  eri^  xsbau  mx/tb^m  suoxoq  b  fix  loiVBrfed  woj  eep 

19d9niBXBq  99'lrfd  9  b  b9bxVOiq  woll  lBXiX9ni-008XV  bflfc  31/008  XV 

.Bisb  •  nilLuze'i  arid  od  wbJ  a'^oxi-.G  pni ylqqs 
-x^q  dn9XBqqs  Id  _  add  bc  vis> ado  pxadns^nxl >1 

.dssd  wold  jn-oq  9lqxdli//n  b  do  ans9m  yd  I>9nxnrx9d9b  sd 

)  ii  o3  noi  ?9X  I jsxdxs. qi-ooa  lv  i  a;  nx  woj  1:  lol 

lj  ooxqx09x  9dd  ei  A'3  dBrf.t  oaxIlBOB?  *£  tdi  >oi  yr  n 


14 


M  ( p i 2  -  p92)  1  w 

— - - -  =  -  +  F  —  (19) 

2y  ZRTL  w/A  k  Ay 

Examination  of  Equation  (19)  indicates  that  a  plot  of  the  left 

_  (27) 

hand  side  versus  w/Ay  ,  termed  a  visco-inertial  flow  plot  , 

is  linear  with  slope  Fn  and  intercept  1/k. 

a 

In  many  instances  application  of  these  graphical 
techniques  yields  anomalous  results  as  evidenced  by  appreci¬ 
able  departures  from  linearity  in  both  the  Klinkenberg  and 
visco-inertial  flow  plots  v  .  In  view  of  the  fact  that  a 
Klinkenberg  extrapolation  is  invalid  if  the  data  are  not  in 
the  viscous  flow  region,  recent  studies  have  revealed  that 
departures  from  linearity  exhibited  by  the  Klinkenberg  plot 

at  high  mean  pressures  may  be  attributed  to  the  inclusion  of 

( 29 ) 

visco-inertial  data  observations  .  Hence,  the  departures 
may  be  avoided,  provided  the  data  points  to  be  used  are  re¬ 
stricted  to  the  viscous  flow  region.  For  flow  in  the  viscous 
region,  examination  of  Equation  (12)  indicates  that  if  per¬ 
meability,  viscosity,  compressibility,  and  temperature  can  be 

assumed  constant  for  a  series  of  data  runs  upon  a  specimen, 

2  2 

a  plot  of  log  (p^  -  p2  )  versus  log  Qq  will  be  linear  with 

slope  of  1.  This  plot  is  termed  the  backpressure  plot.  A 
departure  from  linearity  at  higher  flow  rates  is  assumed  to 
represent  the  existence  of  significant  inertial  effects. 
Consequently,  in  order  to  delineate  the  viscous  and  visco- 
inertial  flow  regions,  Dranchuk  and  Sadiq^2^ 


have  suggested 


. 
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that  data  be  first  plotted  on  a  backpressure  plot.  Data  on 
the  straight  line  portion  of  the  curve  may  be  taken  to  lie 
in  the  viscous  flow  region. 

For  data  taken  in  the  visco-inertial  flow  region, 
downwarping  of  the  data  plot  at  low  Reynolds  Numbers  and 
corresponding  low  mean  pressures  has  been  attributed  to  mole¬ 
cular  streaming.  However,  no  attempt  has  been  made  to  correct 
for  this  effect.  The  approach  in  most  cases  has  been  to  ignore 
the  effect  or  extrapolate  supposedly  unaffected  data. 

Examining  the  Implications  of  Slippage  on  the  Interpretation 

of  Visco-inertial  Flow  Data 

(31) 

Cornell  attributed  the  downwarping  on  the  visco- 

inertial  data  plot  to  molecular  streaming  but  chose  to  avoid 
the  effects  by  extrapolating  unaffected  data.  No  rigid 
criterion  for  delineation  of  the  region  of  linearity  has  been 
established  for  the  visco-inertial  flow  plot.  This  makes  one 
wonder  how  unaffected  data  may  be  chosen.  It  is  equally  likely 
that  some,  all  or  none  of  the  data  are  affected.  Even  though 
slippage  may  contribute  to  the  phenomenon  observed,  the  basic 
theory  of  visco-inertial  flow  may  not  adequately  describe  the 
phenomenon  occurring.  With  this  in  mind,  several  sets  of  the 
data  presented  by  Cornell  were  examined. 
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Results  are  presented  for  two  samples  which  exhibited 
peculiar  characteristic  curves.  Data  observations  had  been 
taken  with  very  similar  gases  so  that  the  Klinkenberg  b-factor 
could  be  taken  as  one  constant  for  each  sample. 

If  the  observed  data  were  consistent  with  the  theory, 
evaluation  of  the  parameters  for  portions  of  the  data  should 
yield  the  same  values  obtained  using  the  entire  data  set. 
Standard  experimental  procedure  is  such  that  the  mean  pressure 
decreases  as  the  flow  rate  decreases.  Thus,  the  effect  of 
slippage  would  become  more  significant  as  the  flow  rate  is 
decreased.  Therefore,  on  a  visco-inertial  flow  plot,  inclusion 
of  data  points  at  successively  lower  values  of  the  abscissa 
w/Ay"  would  increase  the  effect  of  slippage  on  the  determina¬ 
tion  of  permeability  and  the  inertial  resistance  coefficient. 
Hence,  the  data  were  processed  by  first  analyzing  the  two 
highest  flow  rate  data  points  and  successively  adding  the 
lower  flow  rate  data  points,  one  at  a  time,  until  all  of  the 
data  had  been  included.  The  coefficients  k  and  F^,  evaluated 
from  a  linear  least  squares  fit  for  each  group  of  data,  have 
been  plotted  at  the  corresponding  smallest  value  of  w/Ay  in 
the  group  analyzed.  Thus,  this  plot  indicates  the  values 
obtained  for  the  parameters  using  all  available  data  taken  at 
values  of  w/Ay"  above  that  used  for  the  abscissa.  The  differ¬ 
ence  between  any  two  successive  points  indicates  the  effect 
which  the  addition  of  one  more  data  point  at  a  lower  flow  rate 


would  have. 
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(32) 

Results  for  sample  20  are  presented  graphically 

in  Figure  (1) .  The  values  presented  by  Cornell  are  a  per¬ 
meability  k  of  1.01  md  and  a  "turbulence  factor"  F_  of  7.25 

a 

x  10"^  ft  ^ .  However,  it  would  appear  that  at  high  values 

of  the  independent  variable  w/Ay"  the  proper  extrapolated 

values  are  a  permeability  k  of  0.65  md  and  a  "turbulence 

factor"  F  of  5.4  x  10^  ft  ^ .  Generally  it  appears  that 
a 

the  results  for  this  sample  exhibit  internal  inconsistencies. 
It  is  not  immediately  evident  whether  slippage  alone  or  dis¬ 
crepancies  in  the  basic  theory  cause  the  deviation. 

(33 ) 

For  core  sample  14  ,  Cornell  presented  values 

of  20.9  md  for  the  permeability  and  8.2  x  10^  ft  for  the 

"turbulence  factor".  Results  obtained  from  a  linear  least 

squares  fit  of  the  data  are  presented  in  Figure  (2) .  It 

appears  that  extrapolation  should  yield  a  permeability  of 

8  —  1 

16.5  md  and  a  "turbulence  factor"  of  7.0  x  10  ft  * . 

Hence,  the  ramifications  of  not  being  able  to 
distinguish  which,  if  any,  data  are  unaffected  by  slippage 
are  clearly  illustrated  by  indicating  the  error  introduced 
each  time  an  additional  data  point  is  included  in  the  evalu¬ 
ation.  There  is  a  definite  trend  as  the  permeability  and 
"turbulence  factor"  presented  by  Cornell  are  consistently 
higher  than  those  thought  to  be  correct. 

The  theory  indicates  that  both  the  Klinkenberg  and 

visco-inertial  graphical  approaches  should  yield  the  perme- 

.  (34 ) 

ability  of  the  sample.  The  approach  of  Sadiq 


was  to  check 
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20.  OF  CORNELL  USING  STANDARD  VISCO-INERTIAL 
FLOW  PLOT  TECHNIQUE 
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FIGURE  2  FLOW  COEFFICIENTS  FOR  OUTCROP 

SANDSTONE  CORE  14.  OF  CORNELL 
USING  STANDARD  VISCO  -  INERTIAL 
FLOW  PLOT  TECHNIQUE 


xlO'8  (f r1) 
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the  visco-inertial  theory  by  comparing  to  Klinkenberg  plot 
results  of  viscous  data.  He  found  that  generally  the  two  in¬ 
dependent  means  yield  permeabilities  which  differ  by  a  small 
amount.  This  difference  may  be  attributed  to  experimental 
error  or  neglected,  but  close  observation  indicates  that  in¬ 
variably  permeability  obtained  from  the  visco-inertial  plot 
is  high.  When  molecular  streaming  is  pronounced  the  differ¬ 
ence  is  quite  significant. 

For  example,  results  obtained  for  sample  20  of 

c  ^  •  (35) 

Sadiq  were 


k  =  0.000127  md 
and  b  =  160.4  psia 

for  nitrogen  from  a  Klinkenberg  plot  of  the  viscous  data  and 

k  =  0.000184  md 
and  =  8.01  x  1016  ft-1 

from  a  visco-inertial  flow  plot.  Having  observed  a  discre¬ 
pancy  of  45  percent  in  the  permeability  value,  it  is  impos¬ 
sible  to  predict  what  the  error  in  the  other  parameters  may 
be . 

In  the  present  case  and  for  the  others  presented 
the  discrepancies  are  great  enough  to  warrant  performing  a 
correction  for  slippage.  The  delineation  of  unaffected  data 
on  a  visco-inertial  flow  plot  has  been  shown  to  be  very  dif¬ 
ficult,  if  not  impossible,  since  it  is  within  the  realm  of 
possibility  that  in  some  cases  all  data  are  affected. 


'  • 
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Correcting  for  slip  should  enable  the  use  of  all  data  taken 
and  ensure  that  the  analysis  of  viscous  data  and  the  analysis 
of  visco-inertial  data  yield  consistent  results. 
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PROPOSED  THEORY 


Modification  of  the  Visco-inertial  Flow  Theory 


The  Forchheimer  equation 
dp 

-  —  =  F  pq  +  F  pq2  (20) 

dx 

has  been  used  to  describe  visco-inertial  flow  under  the  con¬ 
ditions  of  no  slip.  Klinkenberg ' s  modification  of  Darcy's 
Law  has  also  been  used  to  adequately  describe  the  slip 
phenomenon  under  strictly  viscous  flow  conditions. 

b 

k  =  k(l  +  -)  (21) 

Cl 

P 

The  contribution  of  slip  may  be  incorporated  in  the  viscous 
term  of  the  basic  Forchheimer  equation.  Klinkenberg ' s  modi¬ 
fication  of  Darcy's  Law,  Equation  (10),  may  be  written  as 


dp 

dx 


yq 

k  ( 1  +  -) 
P 


(22) 


From  a  comparison  to  Equation  (9) ,  the  truncated  form  of 
the  Forchheimer  equation,  it  becomes  obvious  that 


1 

k 


k (1  +  -) 
P 


(23) 


Forchheimer ' s  equation  written  in  terms  of  the  absolute  per¬ 


meability  k  is 


dp 

dx 


yq 


k (1  +  -) 
P 


+  FBpq' 


(24) 
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Multiplying  through  by  p/y,  which  is  equal  to 
pM/yZRT  according  to  the  modified  gas  law,  and  defining  the 
constant  mass  flow  rate  as 


w 

pq  =  -  (25) 

A 

results  in 


-M 

- w  pdp 

pZRT  | 


1  fbw 

( - T—  +  — ^-)dx 

k(l  +  |)  Ay 


(26) 


In  order  to  facilitate  the  solution  of  this  equa¬ 
tion  the  apparent  permeability  may  be  taken  as  a  function 
of  the  arithmetic  mean  pressure.  Essentially  this  means 

replacing  p  in  the  viscous  term  by  p.  Substituting  pdp  = 

2 

dp  /2  and  taking  integrals, 

2  p 

-M  P2  dp  L  1  F  w 

-  /  -  =  /  ( - r—  +  — )dx  (27) 

2R  ^  2  yZT  o  k(l  +  =)  yA 

A  p  p 


When  constant  values  of  the  viscosity  and  compres¬ 
sibility,  both  evaluated  at  the  arithmetic  mean  pressure  and 


temperature,  are  used  the  integrated  form  becomes 

1 


M(Pi2  -  Po2) 


2y  ZRTL  ^ 


k  ( 1  +  ±4 
P 


fbw 

a7 


(28) 


At  this  point  it  is  suitable  to  restate  the  Rey¬ 
nolds  Number,  on  the  basis  of  this  development,  as 


N 


Re 


k  FBW 

a7 


b 

(1  +  =) 
P 


(29) 
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It  is  also  appropriate  to  restate  the  integrated 
forms  of  the  viscous  flow  equations  in  the  light  of  Klinken- 
berg ' s  modification  of  Darcy's  Law.  Equations  (11)  and  (12) 
become 

M(px2  -  p,2)  1  1 

—  _  '  w -  =  —  =  - B-  <30) 

2\x  ZRTL  k  k  (1  +  =) 

A  a  p 


A  ToZo(Pl2  -  P22) 


2uZT  P  L  Q 
o  o 


1 

k 


k(!  +  |) 


(31) 


According  to  Equation  (31)  it  is  interesting  to  note  that  a 

b  2  2 

plot  of  log  ( (1  +  =)  (p,  -  p~  ) )  versus  log  Q  ,  which  may  be 

p  J-  ^  o 

termed  a  modified  backpressure  plot,  should  be  linear  of 
slope  1. 


These  relationships  have  been  derived  under  the 
premises  of  a  low  overall  pressure  drop  as  the  questionable 
assumption  that  p  is  not  a  function  of  p  has  been  made. 

To  check  the  ramifications  of  this  assumption  the  theory  was 
rederived  to  account  for  large  pressure  drops.  This  resulted 
in  the  modified  form,  derived  in  detail  in  Appendix  A 


M(p/  ~  P9  ) 


2 u ZRTL  ^ 
A 


Y 


2ab' 


2  2 

k(p,  -p~  ) 


( ab+p, y )  by 

In  - - -  +  — 

(ab+p2y)  kp 


(32) 


where 


and 


fbw 

a  =  — — 

Ay 

1 

y  =  —  +  a 

k 


25 


Methods  for  Solving  the  Modified  Equations 

In  normal  laboratory  analysis,  when  access  to  a 
digital  computer  is  not  available,  it  is  convenient  to  have 
a  graphical  means  of  solution.  Then  a  linear  fit  may  be 
estimated  by  observation  if  necessary. 

Advantage  may  be  taken  of  the  linearity  of  Equation 

w  b 

(28) .  From  a  plot  of  the  left  hand  side  versus  — =  (1  +  =) , 

Ay  p 

referred  to  as  a  modified  visco-inertial  flow  plot,  the  ab¬ 
solute  permeability  is  the  reciprocal  of  the  intercept  and 
the  inertial  resistance  coefficient  is  the  slope.  Use  of 
this  technique  requires  the  prior  knowledge  of  b.  The  modi¬ 
fied  graphical  technique  may  now  be  applied  as  follows: 

1)  The  data  are  plotted  on  a  conventional  backpressure  plot. 

2)  The  viscous  and  visco-inertial  flow  regions  are  ten¬ 
tatively  delineated  by  conservatively  picking  the  lowest 
flow  rate  at  which  departure  from  linearity  appears  to 
occur . 

3)  A  Klinkenberg  plot  is  performed  for  the  data  in  the  vis¬ 
cous  region,  and  tentative  values  of  b  and  k  are  deter¬ 
mined  from  a  linear  least  squares  fit  of  the  data. 

4)  A  modified  backpressure  plot  is  constructed  and  the 
regions  of  viscous  and  visco-inertial  flow  are  again 
delineated.  If  the  point  of  departure  agrees  with  that 
established  in  2)  above,  the  tentative  values  as  deter¬ 
mined  in  3)  are  assumed  correct.  If  not,  the  procedure 
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is  repeated  until  agreement  is  achieved. 

5)  The  final  value  of  b  is  used  to  construct  a  modified 
visco-inertial  plot  for  all  data  except  those  which 
have  been  found  to  be  in  the  viscous  region.  Now  k 
and  F  are  determined  from  a  linear  least  squares  fit 

.D 

of  the  data. 

Graphical  techniques  generally  employed  have  short¬ 
comings  since  the  method  is  time-consuming  and  the  data  must 
be  split  requiring  a  large  number  of  observations  in  both 
the  viscous  and  visco-inertial  flow  regions.  It  would  be 
highly  desirable  to  use  all  of  the  available  data  in  evalu¬ 
ating  each  parameter.  This  would  add  considerably  more 
weight  to  the  values  obtained. 

The  modified  form  of  the  Forchheimer  equation  de¬ 
rived  to  describe  flow  under  relatively  high  pressure  dif¬ 
ferentials  is  not  readily  amenable  to  graphical  treatment. 

Provided  access  to  a  digital  computer  is  avail¬ 
able  a  numerical  procedure  may  be  developed  to  solve  for 
the  best  values  of  the  parameters.  From  an  examination  of 
the  equations  presented  it  is  obvious  that  in  general  the 
relationship  may  be  represented  by 

M(pl2  "  p2 2 } 


2y  ZRTL  £ 
A 


fen (k,b,FB) 


(33) 
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The  functional  relationship  is  non-linear  with  respect  to 
the  three  constants  k,  b  and  Ffi.  Now  it  is  necessary  to  de¬ 
termine  values  of  the  parameters  k,  b  and  FD  to  make  the 
assumed  functional  relationship  optimal  for  experimentally 

M(P]_2  -  P22)  w  _ 

determined  values  of  - -  ■  -  -  ,  — -  ,  p  and  pn  -  p 

2y  ZRTL  ^  Ap7  12 

for  each  of  a  series  of  data  points.  The  criterion  of  op¬ 
timization  which  may  be  used  for  the  estimation  of  the  best 
fit  parameters  is  minimization  of  the  residual  sum  of  squares 
of  the  left  hand  side. 

Since  this  type  of  solution  hinges  upon  linearity, 

quasilinearization,  a  combination  of  linear  approximation 

techniques  with  the  capabilities  of  the  digital  computer, 

(36) 

may  be  applied.  The  Newton-Raphson  approximation  scheme v  , 
by  analogy  with  the  sequence  of  approximations  to  the  zeros 
of  a  function,  may  be  applied  to  present  an  expansion  about 
a  set  of  parameter  values.  This  results  in  the  iterative 
solution  of  a  set  of  linear  equations.  Thus,  given  initial 
estimates  one  may  apply  linear  least  squares  to  determine 
corrections  to  the  guesses.  Convergence  is  not  always 
achieved  and  is  related  to  the  selection  of  initial  guesses. 
If  the  initial  parameter  estimates  are  good,  convergence 
will  occur  rapidly. 

Application  of  this  technique,  which  is  described 
in  detail  in  Appendix  B,  requires  evaluation  of  several 
derivatives  of  the  functional  relationship  with  respect  to 
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the  parameters.  In  Appendix  C  the  necessary  relationships 
are  developed  for  the  two  equations  presented  previously. 

Estimation  of  the  Parameters 

In  order  to  assure  convergence  in  the  interpreta¬ 
tion  of  data  by  the  numerical  technique  proposed  it  is 
extremely  important  to  have  good  initial  guesses  of  the 
unknown  parameters.  In  a  re-evaluation  of  the  results  of 
previous  studies,  the  values  presented  in  the  original  work 
leave  no  guess  work  for  good  initial  guesses.  However,  in 
the  case  of  recent  tests,  where  only  basic  physical  data 
are  available,  a  simple  technique  must  be  determined  for 
making  good  estimations  of  the  parameters. 

Permeability  may  be  approximated  from  a  single 
point  test.  Any  single  point  could  do  but  it  is  suggested 
the  lowest  flow  rate  observation  be  used. 

The  Klinkenberg  b-factor  for  air  has  been  cor- 

( 37 ) 

related  with  absolute  permeability  by  Heid  et  al  as 

b  =  0.777  k  °  *  (34) 

for  b  in  international  atmospheres  and  k  in  millidarcies . 
It  is  suggested  that  b,  calculated  from  this  equation  using 
the  single  point  permeability,  be  used  as  a  starting  value. 
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(38) 

As  reported  by  Sadiq  many  attempts  have  been 

made  to  predict  the  inertial  resistance  coefficient.  From  a 
dimensional  analysis  he  proposed 


1.617  x  105 

FB  =  .  .5  .5.07 

k  <j> 


(35) 


where  Fg  is  in  ft  ^  and  k  in  md,  as  a  general  correlation  of 
Fd  with  the  absolute  permeability  and  the  effective  porosity. 
It  is  therefore  suggested  that  this  relationship  be  used  to 
obtain  a  good  initial  guess  of  Fg  using  the  single  point 
permeability  evaluation. 
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EXPERIMENTAL  APPARATUS 

A  gas  expansion  porosimeter,  calibrated  with  solid 
blanks  of  known  dimensions,  was  used  to  measure  the  gross 
grain  volume  of  the  cores  prior  to  mounting. 

Accurate  determination  of  the  characteristic  flow 
parameters  of  the  samples  tested  required  precise  measurement 
and  control  of  temperature,  pressure  and  volumetric  flow  rate 
over  an  extensive  range  of  values  for  each  data  set. 

The  temperature  of  the  flowing  fluid  immediately 
upstream  and  downstream  of  the  endface  flow  taps  was  measured 
with  calibrated  iron  constantan  thermocouples,  inserted  in 
the  flow  line.  Control  of  the  mean  temperature  was  effected 
by  passing  the  fluid  through  a  copper  coil,  immersed  in  a 
constant  temperature  glycol  bath,  before  entering  the  sample. 

Water  and  mercury  manometers  were  used  to  obtain 
pressure  readings  up  to  30  psig,  whereas  higher  pressures 
were  measured  with  a  variety  of  calibrated  bourdon-tube 
pressure  gauges.  Control  of  the  upstream  pressure  was  ef¬ 
fected  by  means  of  pressure  regulators  in  series  between  the 
gas  supply  cylinders  and  the  sample  holder.  Downstream  pres¬ 
sure  could  be  controlled  by  a  backpressure  regulator  in  the 
flow  line. 

Several  calibrated  volumetric  flow  measuring  devices 
were  used  so  that  an  extensive  range  could  be  covered  ade¬ 
quately.  Flow  rates  below  15  cc/sec  were  measured  by  timing 


■ 
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the  displacement  of  soap-film  in  a  burette.  Rates  between 
15  cc/sec  and  180  cc/sec  were  measured  with  a  wet  test  meter 
and  higher  rates  to  1650  cc/sec  were  measured  using  a  dia¬ 
phragm  meter. 

The  coreholders  used  were  of  a  rubber  sleeve  con¬ 
struction  and  a  resin  mounted  type.  For  the  sleeve  type  of 
core  holder  illustrated  in  Figure  (3)  a  1/8-inch  neoprene 
sleeve  was  used  satisfactorily  with  nitrogen  as  the  pressurizing 
fluid,  effecting  a  seal  around  the  core  and  maintaining  a 
longitudinal  compression  on  the  core.  The  core  holder  was 
modified  by  penetrating  the  rubber  sleeve  with  1/8  inch 
stainless  steel  Autoclave  tubing  to  enable  the  measurement  of 
the  pressure  distribution  along  the  length  of  the  sample. 

(39) 

The  resin  mounted  core  holder  has  been  described  by  Mackett 

The  inlet  and  outlet  flow  lines  and  all  pressure  taps 
consisted  of  1/8-inch  stainless  steel  Autoclave  tubing.  Tygon 
tubing  was  used  to  connect  the  flow  measuring  devices.  The 
entire  apparatus  was  enclosed  in  a  thermostatically  controlled 
constant  temperature  cabinet. 
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EXPERIMENTAL  PROCEDURE 

The  core  samples  were  cleaned  by  refluxing  toluene 
vapors  at  240°F  for  20  hours  and  drying  at  220°F  for  24 
hours . 

The  bulk  volume  of  each  of  the  cores  was  evaluated 
from  extensive  measurements  of  length  and  diameter.  The 
effective  grain  volume  of  each  sample  was  determined  using  a 
calibrated  gas  expansion  porosimeter.  From  these  measure¬ 
ments  the  effective  porosity  was  evaluated. 

The  cores  were  subsequently  mounted  in  the  sleeve 
type  holder  or  mounted  in  resin  in  a  similar  fashion  to  that 
proposed  by  Mackett,  except  that  the  filler  in  the  Hysol  4160 
resin  was  excluded. 

Following  the  insertion  of  the  core  into  the  core 
holder  and  the  connection  of  the  pressure  taps  and  flow  lines, 
the  apparatus  was  pressure  tested  for  leaks  using  a  soap 
solution  type  detector. 

Subsequently,  the  desired  flow  tests  were  run,  with 
data  points  chosen  in  a  random  order  as  a  check  for  any  in¬ 
ternal  inconsistencies.  For  all  tests  on  the  sleeve  type 
core  holder  the  pressurizing  fluid  pressure  was  maintained  100 
psi  above  the  inlet  pressure.  The  temperature  of  the  cabinet 
was  observed  periodically  to  ensure  that  it  remained  constant 
during  a  series  of  runs  and  the  temperature  of  the  glycol  bath 
was  adjusted  so  that  the  arithmetic  mean  temperature  of  the 
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flowing  fluid,  nitrogen,  was  constant. 

For  each  run,  after  adjustment  of  the  pressure 
regulating  equipment,  when  observation  of  the  inlet  and  out¬ 
let  pressures  and  volumetric  flow  rate  indicated  that 
stabilized  conditions  prevailed,  values  were  recorded  of 
1)  barometric  pressure  and  ambient  temperature,  2)  inlet  and 
outlet  thermocouple  readings,  3)  inlet,  outlet  and  inter¬ 
mediate  pressure  readings,  and  4)  volumetric  flow  rate  reading 
with  the  corresponding  pressure  and  temperature  readings  for 
the  particular  measuring  device. 

For  those  series  of  runs  in  which  a  resident  uni¬ 
form  liquid  saturation  was  desired  2-octanol,  a  chemically 
inactive  fluid  possessing  a  low  vapor  pressure,  was  used  as 
the  saturating  fluid.  In  order  that  the  proper  proportions 
of  this  fluid  be  deposited  within  the  core  it  was  found  neces¬ 
sary  to  mix  with  a  non-reacting  solvent  of  high  vapor  pressure, 
such  as  pentane,  which  could  easily  be  removed  from  the  core. 

Hence,  the  entire  procedure  consisted  of  first 
evacuating  the  core  sample  in  situ.  Then  a  deaerated  solu¬ 
tion  of  the  fluids,  containing  the  saturating  fluid  in  the 
desired  proportion,  was  injected  into  the  sample  with  a  Ruska 
proportioning  pump  to  a  pressure  higher  than  that  encountered 
under  flow  conditions,  noting  in  each  case  the  volume  required 
to  fill  the  holder  and  sample.  Saturations  required  were 
generally  below  the  minimum  saturation  for  liquid  flow,  hence, 
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an  immobile  saturation  of  mixture  was  established  by  flowing 
in  both  directions  at  a  high  pressure  and  a  considerably 
higher  rate  than  those  expected  under  test  conditions.  Then 
the  solvent  was  removed  by  flowing  at  low  pressures  and  under 
vacuum  for  sufficient  periods  of  time.  Between  series  of 
runs  the  sample  was  removed  from  the  coreholder  and  weighed 
to  determine  the  resident  liquid  saturation. 
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TREATMENT  OF  DATA 


A  general  computer  program  was  written  to  process 
the  recorded  experimental  gas  flow  data.  Observed  pressure 
readings  and  volumetric  flow  rate  readings  were  corrected 
and  reduced  to  standard  conditions  by  applying  the  appropri¬ 
ate  conversion  and  calibration  factors.  Volumetric  flow  rate 
was  recorded  at  the  reference  conditions  of  1  international 
atmosphere  and  520°R. 

In  determining  fluid  properties  the  values  of  vis¬ 
cosity  and  compressibility  were  evaluated  at  the  arithmetic 
mean  temperature  and  pressure.  The  relationship  used  for 
the  viscosity  of  nitrogen  was  that  suggested  by  Kestin  and 


M  =  1778  x  10"7(1  +  8.958  x  10_4  (p-1) 

+  6.120  x  10"7  (p-1)2  +  3.997  x  10~8(p-l)3) 
+  4.55  x  10~7  (T  -  25)  (36) 


for  viscosity  y  in  poise,  pressure  p  in  international  atmos¬ 
pheres  and  temperature  T  in  °C . 

Compressibility  was  evaluated  by  four-point  Lagran- 

(41) 

gian  interpolation  of  curve  fitted  data  of  Hilsenrath  et  al 
at  temperatures  of  280  ,  290  ,  300  ,  and  310°K. 

The  expressions  necessary  for  performing  graphical 
plots  and  mathematical  fits  of  the  processed  data  were  recorded 
on  IBM  cards  to  facilitate  their  use  in  further  processing. 
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Programs  were  also  written  to  perform  least  squares 
fit  of  viscous  data  on  a  Klinkenberg  plot,  least  squares  fit 
of  Forchheimer  equation  visco-inertial  flow  plot  and  modified 
visco-inertial  flow  plot,  quasilinearization  of  the  modified 
visco-inertial  flow  equation  for  low  pressure  gradients  and 
the  modified  visco-inertial  flow  equation  for  high  pressure 
gradients.  The  quasilinearization  technique  was  programmed 
for  an  IBM  7040  using  the  method  of  determinants  and  Cramers 
Rule  for  the  simultaneous  solution  of  the  linear  equations. 

A  sample  of  the  program  and  related  nomenclature  appears  in 
Appendix  C. 
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RESULTS 


Reanalysis  of  the  Data  of  Previous  Investigators 

In  order  to  apply  the  modified  visco-inertial  flow 
graphical  technique  it  is  necessary  to  have  a  value  of  the 
Klinkenberg  b-factor.  For  sample  20  of  Sadiq  a  Klinkenberg 
extrapolation  had  been  performed  on  the  viscous  data  and 
b  had  been  evaluated  as  160.4  psi  for  nitrogen.  Therefore, 
the  modified  technique  was  applied  to  the  visco-inertial 
flow  data  of  sample  20^2^. 

The  numerical  technique  was  applied  to  both  modi¬ 
fied  visco-inertial  flow  equations,  using  the  original  de- 

(43) 

terminations  for  initial  guesses  .  A  summary  of  the  re¬ 
sults  obtained  is  presented  in  Table  1. 

Table  1 


A  Summary  of 

Results  for 

Core  20  of 

Sadiq 

Method 

k 

(md) 

b 

(psi) 

F 

B_ 

(ft 

Klinkenberg  Plot 

0.  000127 

160.4 

Visco-inertial 

Flow  Plot 

0.000184 

8.01  x 

1016 

Modified  Visco- 
inertial  Flow  Plot 

0.  000124 

(160.4) 

2.58  x 

1016 

Quasi linearization 
of  Eqn.  (28) 

0.  000125 

172.3 

3.18  x 

1016 

Quasi linearization 
of  Eqn.  (32) 

0.  000126 

171.3 

3.22  x 

1016 
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Since  sufficient  viscous  data  were  not  available 
for  the  two  samples  of  Cornell,  quasilinearization  was  per¬ 
formed  on  both  modified  equations.  The  values  presented  by 
Cornell  were  used  as  initial  guesses  and  b  was  estimated 
using  the  correlation  previously  proposed.  A  summary  of  re¬ 
sults  obtained  for  samples  20  and  14  is  presented  in  Table 
2. 


Table  2 


Quasilinearization  Results  for 
Samples  20  and  14  of  Cornell 


Method 

k 

(md) 

b 

(psi) 

fb 

(ft 

Sample  20 

Cornell 

1.01 

7.25  x 

10 

Estimation 

0.65 

5.4  x 

10 

Quasi linearization 
of  Eqn.  (28) 

0  .595 

67.5 

5.66  x 

10 

Quasi linearization 
of  Eqn .  ( 32 ) 

0  .614 

65.0 

5.75  x 

10 

Sample  14 

Cornell 

20.9 

8.2  x 

10 

Estimation 

16.5 

7.0  x 

10 

Quasi linearization 
of  Eqn.  (28) 

18.28 

11.97 

7.79  x 

10 

Quasi linearization 
of  Eqn.  (32) 

18.35 

11.89 

7.8  x 

10 

10 

10 

10 

10 
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Present  Experimental  Results 

Sufficient  data  were  taken  for  several  core  samples 
so  that  a  careful  evaluation  of  the  basic  theory  and  the  slip 
correction  for  flow  both  in  a  dry  core  sample  and  in  the  pre¬ 
sence  of  an  immobile  liquid  phase  could  be  made.  The  physical 
properties  of  these  samples  are  listed  in  Table  3. 

Table  3 

Physical  Properties  of  the  Samples 


Sample 

1 

2 

3 

Description 

Berea 

Sandstone 

Gilwood 

Sandstone 

Swan  Hills 
Limestone 

Diameter 

(cm) 

8.93 

8.82 

8.66 

Length 

(cm) 

21.8 

18.3 

19.9 

Area 

,  2 » 

(cm  ) 

62.7 

61.2 

58.9 

Pore  Volume 

(cm  ) 

293. 

276. 

165. 

Porosity 

0.214 

0.246 

0.141 

A 

low  permeability  berea 

sandstone , 

sample  1,  was 

chosen  for  a 

series 

of  runs  because 

of  the  probable  homogeneity 

of  the  sample.  The  basic  flow  data  appear  in  Table  10,  Appen¬ 
dix  D  and  are  presented  graphically  on  a  standard  backpressure 
plot  in  Figure  (4).  As  expected,  the  viscous  data  do  not  lie 
on  a  45°  straight  line.  The  departure  from  linearity  and 


41 


"QO"or  Qo  (cc/sec) 


FIGURE  4 


BACKPRESSURE  PLOT  FOR  CORE  1 
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hence  the  transition  between  viscous  and  visco-inertial  flow 
appears  to  occur  at  a  rate  of  4  cc/sec.  A  least  squares  fit 
of  the  tentative  viscous  data  observations,  presented  on  a 
Klinkenberg  plot  in  Figure  (5) ,  resulted  in  values  of 

k  =  2.29  md 

and  b  =  0.610  atm.  for  nitrogen. 

As  a  check  for  internal  inconsistencies  k  and  b  were 
calculated  using  the  six  observations  at  the  lowest  mean  pres¬ 
sures  and  thereafter  for  groups  of  observations  formed  by 
successively  adding  one  point  at  the  next  higher  mean  pres¬ 
sure.  The  results  were  plotted  in  Figure  (6)  with  the  ab¬ 
scissa  being  the  lowest  reciprocal  mean  pressure  in  each 
group.  The  data  appear  to  be  quite  consistent  and  the  re¬ 
sults  presented  for  the  tentative  viscous  data  appear  to  be 
quite  representative. 

A  standard  visco-inertial  flow  plot  was  performed 

for  all  data.  It  appeared  that  the  data,  plotted  in  Figure 

(7),  could  be  fitted  by  a  straight  line  for  w/Ay  greater 

2 

than  0.3  gm/sec  cm  cp .  These  data  could  be  assumed  un¬ 
affected  by  slippage.  Extrapolation  of  these  data  resulted 
in  values  of 

k  =  2 . 2  0  md 

and  F_.  =  143.  atm  sec^/gm  =  4.42  x  lO^ft  ^ 

a 

A  check  was  made  using  all  nonviscous  data,  a  linear  fit 
yielding  values  of 
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FIGURE  5  KLINKENBERG  PERMEABILITY  PLOT 

FOR  CORE  1. 
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FIGURE  6  KLINKENBERG  PERMEABILITY  PLOT  COEFFICIENTS  FOR 

CORE  1. 
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FIGURE  7  STANDARD  VISCO- INERTIAL  FLOW  PLOT  FOR  CORE 
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k  =  2.41  md 

and  Fg  =  225.  atm  sec^/gm  =  6.95  x  10^  ft  ^ 

The  results  of  an  analysis,  identical  to  that  per¬ 
formed  upon  Cornell's  data,  are  presented  in  Figure  (8).  As 
noted  previously,  the  results  seem  to  be  quite  consistent  for 
w/Ay  values  larger  than  0.3.  Nevertheless  there  is  a  definite 
trend  as  the  increasing  effect  of  slippage  causes  the  values 
of  k  and  F  to  increase  continuously  with  decreasing  w/Ay . 

Using  the  b  value  of  .610  atm,  determined  from  the 
Klinkenberg  plot,  a  modified  visco-inertial  flow  plot  was 
made  and  is  presented  in  Figure  (9).  A  linear  least  squares 
fit  gave  values  of 

k  =  2.26  md 

and  F-  =  184.  atm  sec^/gm  =  5.69  x  10^  ft  ^ 

b 

As  previously,  data  were  processed  beginning  at  the 

highest  flow  rates.  The  results  are  plotted  in  Figure  (10) 

with  the  abscissa  (1  +  —  )w/AyT.  The  modified  analysis  tech- 

F 

nique  does  not  seem  to  correlate  the  data  very  well  although 
the  results  appear  to  be  more  consistent  than  those  obtained 
using  the  standard  technique. 

From  the  proposal  made  previously  the  initial  guesses 
for  the  coefficients  of  this  sample  were  evaluated  as 

k  =  3 . 69  md 

b  =  0.467  atm 

F_.  =  6.66  atm  sec^/gm  =  2.06  x  10^  ft 


and 


. 


k  ( md ) 
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FIGURE  8  FLOW  COEFFICIENTS  FOR  CORE  1.  USING 

STANDARD  VISCO- INERTIAL  FLOW  PLOT 
TECHNIQUE 
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FIGURE  9  MODIFIED  VISCO- INERTIAL  FLOW  PLOT  FOR  CORE  1. 
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FIGURE  10  MODIFIED  VISCO  -  INERTIAL  FLOW  EQUATION  COEFFICIENTS 

FOR  CORE  1. 
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Quasilinearization  was  performed  on  both  Equations  (28)  and 

-4 

(32) .  With  a  convergence  criterion  e  =  10  ,6  iterations 

were  required  in  each  case  to  converge  to  similar  results, 
identical  to  three  significant  figures, 

k  =  2.24  md 

b  =  .670  atm 

2  9 

and  Fg  =  177.  atm  sec  /gm  =  5.47  x  10  ft 

The  results  predicted  by  using  these  coefficients  are  repre¬ 
sented  by  the  dashed  line  on  the  visco-inertial  flow  plot, 
Figure  (7) . 

In  the  graphical  analysis  it  is  necessary  to  assume 

that  b  is  a  constant.  Using  the  numerical  technique  it  is 

possible  to  examine  the  constancy  of  b  as  well  as  k  and  F_ . 

a 

A  repetitive  analysis  was  performed  in  solving  Equation  (32) 
numerically,  by  successively  using  more  data  points  in  the 
direction  of  lower  w/Ay  .  Results  plotted  in  Figure  (11) 
against  the  abscissa  w/Ay  indicate  that  there  is  no  consis¬ 
tency  until  all  of  the  data  has  been  used.  A  check  made 
using  series  of  four  representative  data  points  indicated 
that  good  results  may  be  obtained  provided  the  data  observa¬ 
tions  fully  cover  the  flow  regions  investigated. 

A  summary  of  the  results  obtained  for  Core  1  is 


listed  in  Table  4. 
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FIGURE  11  RESULTS  FROM  QUASILINEARIZATION  OF  EQUATION  (32) 

FOR  CORE  1. 
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Table  4 

A  Summary  of  Results  for  Core  1 


Method 

k 

(md) 

b 

( atm . ) 

atm  sec^ 

gm 

fb 

ft-1 

Klinkenberg  Plot 

2.29 

0.610 

Visco-Inertial  Flow 

Plot 

2.20 

143. 

4.42  x 

io9 

Modified  Visco- 

Inertial  Flow  Plot 

2.26 

(.610) 

184. 

5.69  x 

109 

Quasilinearization 

of  Eqn.  (38) 

2.24 

0.670 

177. 

5.47  x 

109 

Quasi linearization 

of  Eqn.  (32) 

2.24 

0.670 

177. 

5.47  x 

109 

Now  using  the  best 

values,  obtained  from  the 

numerical  analysis  using  all 

data,  a  check  on  the  cut-off 

point  indicated  a  Reynolds  Number  of  .00199. 

In  order  to  check  for  homogeneity  of  the  sample  and 
the  presence  of  mechanical  end  effects,  a  pressure  drop  pro¬ 
file  was  taken  during  several  of  the  runs  on  this  sample. 
Graphical  results  of  the  data  in  Tables  11  to  14  are  pre¬ 
sented  in  Figure  (12) .  The  results  obtained  using  b  =  0.610 
atm  are  presented  in  Table  5. 
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PRESSURE  DROP  PROFILE  FOR  CORE 
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Table  5 

Permeability  and  Inertial  Resistance  Coefficient 

Distribution  for  Core  1 


2 

(md)  (atm  sec  /gm) 


Section 

1 

2.15 

340  . 

Section 

2 

2.10 

200  . 

Section 

3 

2.44 

255  . 

Section 

4 

2.50 

290  . 

An  analysis,  similar  to  that  performed  for  sample 
1,  was  carried  out  upon  the  data  for  the  Gilwood  sandstone 
sample  number  2.  Experimental  gas  flow  data  are  presented 
in  Table  15,  Appendix  D.  From  a  backpressure  plot  of  these 
data,  Figure  (13),  a  straight  line  may  be  fitted  to  the  data 
obtained  at  flow  rates  less  than  40  cc/sec.  under  reference 
conditions.  A  fit  of  the  tentative  viscous  data  on  a  Klin- 
kenberg  plot,  Figure  (14),  indicated 

k  =  3  6 . 5  md 

and  b  =  0.401  atm.  for  nitrogen 

An  analysis  made  by  successively  including  more  data  points 
in  the  evaluation,  presented  graphically  in  Figure  (15), 
illustrates  that  the  data  are  quite  consistent  internally 
and  that  the  tentative  results  are  representative. 
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The  apparently  straight  line  portion  of  the  standard 
visco-inertial  flow  plot  in  Figure  (16)  occurs  at  values  of 
w/Ay  greater  than  0.7.  Results  obtained  from  a  linear  least 


squares  fit  were 

k 

=  35.8  md 

and  Fd 

a 

=  4.80  atm  sec2/gm  =  1.48  x  108  ft  1 

Using  the  modified  visco-inertial  technique  illus¬ 
trated  in  Figure  (17),  values  obtained  using  the  nonviscous 


data  were 

k 

=  36.4  md 

and  F 

Jd 

=  5.78  atm  sec2/gm  =  1.79  x  108  ft  ^ 

Results  of  analyses,  similar  to  those  performed  upon  the 
data  of  sample  1,  are  presented  graphically  in  Figures  (18) 
and  (19)  . 

For  the  quasilinearization  schemes,  initial  guesses 
were  obtained  using  the  correlations  previously  presented 


k 

=  51.5  md 

b 

=  0.166  atm 

and  F 

J3 

=  0.874  atm  sec2/gm  =  2.70  x  10^  ft  ^ 

-4 

For  a  convergence  criterion  e  of  10  ,  solution 

of  Equations ( 28 )  and  (32)  both  yielded  the  results 


k 

=  36.2  md 

b 

=  0.423  atm  for  nitrogen 

and  Fg 

=  5.62  atm  sec2/gm  =  1.74  x  108  ft  ^ 

and 
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FIGURE  16  STANDARD  VISCO-INERTIAL  FLOW  PLOT  FOR  CORE  2 
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exact  to  three  significant  figures,  with  6  iterations  re¬ 
quired  for  the  solution  of  Equation  (28)  and  9  iterations  for 
the  solution  of  Equation  (32). 

A  summary  of  the  results  obtained  for  core  2  is  pre¬ 
sented  in  Table  6 . 


Table  6 

A  Summary  of  Results  for  Core  2 


Method 

k 

(md) 

b 

(atm) 

,atm  sec^, 
gm 

(ft" 

h 

Klinkenberg  Plot 

36.5 

0.401 

Vi sco- Inertial 

Flow  Plot 

35.8 

4.80 

1.48  x 

108 

Modified  Visco- 
Inertial  Flow  Plot 

36.4 

(.401) 

5.78 

1.79  x 

108 

Quasi linearization 
of  Eqn.  (28) 

36.2 

0.423 

5.62 

1.74  x 

108 

Quasi linearization 
of  Eqn.  (32) 

36.2 

0.423 

5.62 

1.74  x 

108 

Using  the  values  obtained  by  quasilinearization  a 
check  of  the  Reynolds  Number  at  the  cut-off  indicated  a 
value  of  0.  012  . 

Since  this  sample  was  a  vertical  core  from  a  produc¬ 
ing  formation,  it  was  found  necessary  to  check  the  degree  of 
homogeneity.  Using  pressure  taps  installed  along  the  length 
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of  the  sample,  several  pressure  profiles  were  measured.  The 
data  in  Tables  16  to  19,  Appendix  D  are  exhibited  graphically 
in  Figure  (20).  Analysis  of  the  data  for  each  section,  using 
the  numerical  techniques,  resulted  in  the  average  values 
listed  in  Table  7. 


Table  7 

Permeability  and  Inertial  Resistance 
Coefficient  Distribution  for  Core  2 


F 

kb  B 

(md)  (atm)  (atm  sec2/gm) 


Section 

1 

13.2 

.917 

9.76 

Section 

2 

92.2 

.565 

2.13 

Section 

3 

254. 

.208 

.894 

Section 

4 

124. 

.452 

1.07 

In  order  to  check  that  the  effect  observed  was  due 
to  heterogeneity  rather  than  an  entrance  effect,  the  direc¬ 
tion  of  flow  was  reversed  and  a  test  was  performed  simulat¬ 
ing  one  of  the  flow  rates.  The  pressure  profiles  are  com¬ 
pared  in  Figure  (21) . 

On  the  premises  that  the  numerical  technique  is 
vastly  superior,  the  applicability  of  the  visco-inertial 
theory  for  flow  in  the  presence  of  a  liquid  phase  was  examined. 
Data  for  the  Swan  Hills  limestone  core  sample  3  were  taken  at 
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67 


four  immobile  liquid  saturation  conditions  as  well  as  for 
the  dry  core. 

Basic  flow  data  for  the  core,  dry  and  at  liquid 
saturations  of  0.0765,  0.124,  0.205,  and  0.253,  are  pre¬ 
sented  in  Tables  20  to  24,  Appendix  D.  These  data  were 
subjected  to  the  numerical  solution  of  Equation  (32)  re¬ 
sulting  in  the  values  presented  in  Table  8.  The  solution 
of  Equation  (28)  yielded  very  similar  results. 

Table  8 

Effective  Permeability,  Klinkenberg  b-Factor  and 
Inertial  Resistance  Coefficient  as  a  Function 
of  Liquid  Saturation  for  Core  3 

Liquid 

Saturation  k  b 

S  (md)  (psia) 


Dry 

13.0 

0.649 

85.5 

2.64 

X 

10 

0.0765 

11.27 

0.608 

101. 

3.12 

X 

10 

0.124 

10.0 

0.501 

125. 

3.86 

X 

10 

0.205 

8.17 

0.611 

148. 

4.57 

X 

10 

0.253 

6.36 

0.574 

228  . 

7.05 

X 

10 

( 


atm  sec' 
gm 


B 


(ft  1) 


The  results  predicted  by  using  the  calculated 
values  of  the  three  parameters  are  compared  with  the  physi¬ 
cal  observations  on  the  visco-inertial  flow  plot.  Figure  (22). 

The  coefficients  are  plotted  as  a  function  of  liquid  saturation  in 


LIQUID  SATURATIONS 
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Figure  (23) .  The  effective  permeability  was  found  to  de¬ 
crease  steadily  with  every  increase  in  the  liquid  saturation. 
Also  the  Klinkenberg  b-factor  was  found  to  be  relatively 
constant.  The  inertial  resistance  coefficient  was  found  to 
increase  with  liquid  saturation  at  an  increasing  rate. 


Modification  of  the  Correlation  of  the  Inertial  Resistance 

Coefficient 


When  the  values  determined  from  Sadiq's  correlation 

for  the  inertial  resistance  coefficient  were  compared  to  the 

experimentally  determined  values  of  a  large  discrepancy 

was  noted.  A  check  on  Sadiq's  approach  revealed  that  his 

proposal  had  been  formed  using  dimensional  analysis  based 

on  semi-theoretical  reasoning,  and  involved  a  least  squares 

(44) 


fit  of  the  data  of  Cornell 


„  . ,  .  (45)  ,  _  ,  .  (46) 

Hamilton  ,  and  Sadiq 


It  cannot  be  assumed  that  the  three  parameters 
effective  porosity,  absolute  permeability  and  the  inertial 
resistance  coefficient  should  completely  characterize  the 
structure  of  a  porous  medium.  Nevertheless,  if  these  factors 
are  assumed  to  be  related  a  least  squares  fit  of  the  basic 
data  should  yield  the  best  predictive  tool. 

The  relationship  which  has  been  shown  to  have  some 
theoretical  basis  is 


c  ,  d 

Fb  =  a  k 


(37) 


where  a,  c  and  d  are  constants. 
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FIGURE  23  EFFECTIVE  PERMEABILITY,  KLINKENBERG 

b-FACTOR  AND  INERTIAL  RESISTANCE 
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This  was  linearized  by  taking  logarithms 


In  Fg  =  In  a  +  c  In  cf>  +  d  In  k 


(38) 


A  least  squares  fit  for  Cornell's  basic  data  yielded 

2.80  x  109 

FB  =  1.09  .1.17 

<p  k 


where  F  is  in  ft  ^  and  k  is  in  millidarcies .  Hamilton's 
a 

data  were  subjected  to  a  similar  analysis  resulting  in 


6.56  x  108 

FB  =  1.56  .1.22 

4>  k 

while  a  fit  of  Sadiq's  data  yielded 

3.40  x  104 

FB  =  6.31  ,0.50 

<J>  k 

A  composite  fit  of  all  the  data  indicated 

7.56  x  108 

FB  =  1.67  .1.12 

<p  k 


(40) 


(41) 


(42) 


A  comparison  of  results  obtained  using  the  proposed 
techniques  with  the  experimental  results  is  given  in  Table  9. 
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Table  9 

A  Comparison  of  Inertial  Resistance 
Coefficient  Estimations 


Experimental 
Sample  Value 


Fb  (ff1 2 3) 

Sadiq ' s 
Correlation 


1  5.47  x  109  2.64  x  108 

2  1.74  x  108  3.2  x  107 

3  2.64  x  109  9.30  x  108 


Present 

Correlation 

3.98  x  109 
1.39  x  108 
1.13  x  109 


DISCUSSION  OF  RESULTS 


For  core  20  of  Sadiq  the  permeability  value  obtained 
by  the  modified  visco-inertial  flow  plot  technique  is  con¬ 
sistent  with  that  obtained  from  the  Klinkenberg  extrapolation 
of  the  viscous  data.  However,  the  value  of  the  inertial  re¬ 
sistance  coefficient  is  approximately  one-third  of  the  value 
obtained  using  the  standard  visco-inertial  flow  plot  technique. 
The  results  obtained  by  quasilinearization  of  the  two  modified 
equations  compare  very  closely.  They  are  quite  consistent 
with  the  results  obtained  from  the  Klinkenberg  and  modified 
visco-inertial  flow  plots. 

One  important  observation  which  may  be  made  about 
the  data  is  that  not  one  data  point  can  be  considered  free 
of  slippage  effects.  Even  at  the  highest  mean  pressure  con¬ 
dition  the  factor  (1  +  b/p)  is  found  to  be  greater  than  1.20. 
The  significance  of  this  observation  is  that  the  highest  flow 
rate  data  have  been  obtained  at  the  extreme  pressure  capacity 
of  the  equipment.  Hence,  the  equipment  would  have  to  be 
modified  especially  for  this  sample  if  unaffected  data  were 
required.  This  indicates  the  superiority  of  schemes  whereby 
slippage,  be  it  especially  significant  or  not,  can  be  ac¬ 
counted  for. 

For  the  two  cases  in  which  the  data  of  Cornell  were 
used  the  lack  of  sufficient  viscous  data  for  a  Klinkenberg 
extrapolation  prevented  the  evaluation  of  the  Klinkenberg 
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b-factor  and  hence  prevented  the  use  of  the  modified  visco- 
inertial  flow  plot  technique.  Since  most  of  the  data  were 
affected  by  slippage,  only  very  few  of  the  points  could  be 
used  to  determine  k  and  F_.  using  the  standard  visco-inertial 
flow  plot.  Hence,  the  numerical  technique  was  the  only 
feasible  means  of  solution.  Application  of  the  numerical 
technique  to  the  modified  forms  of  the  visco-inertial  flow 
equation  gave  results  which  were  very  consistent  with  those 
estimated  by  extrapolation  of  the  visco-inertial  flow  plot. 

The  large  deviations  found  when  a  comparison  with 
Cornell's  values  was  performed  casts  doubt  upon  the  reli¬ 
ability  of  results  presented  by  previous  investigators  who 
have  used  the  standard  technique. 

In  the  present  work,  since  very  many  data  observa¬ 
tions  had  been  taken  for  both  core  samples  1  and  2,  it  was 
possible  to  distinguish  a  region  of  linearity  on  a  back¬ 
pressure  plot  quite  readily.  As  expected,  the  viscous  data 
did  not  lie  on  a  45°  slope  in  both  cases.  Generally  the 
tentative  viscous  results  were  found  to  be  quite  representa¬ 
tive,  as  indicated  by  Figures  (6)  and  (15) .  These  plots  also 
indicate  that  the  addition  of  several  data  points  from  the 
visco-inertial  flow  region  may  cause  serious  errors  in  the 
values  obtained  for  the  absolute  permeability  and  the  Klinken- 
berg  b-factor. 

In  both  cases  a  region  of  linearity  could  also  be 
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distinguished  on  a  standard  visco-inertial  flow  plot.  However, 
Figures (8)  and  (18)  clearly  illustrate  that  within  the  region 
of  linearity  the  values  of  the  parameters  k  and  F  consis- 
tently  increase  with  decreasing  w/A jT.  The  plots  also  indicate 
that  the  use  of  several  affected  data  points  can  cause  serious 
errors  in  the  calculated  values  of  k  and  F  .  Generally  the 
values  of  k  did  not  correspond  very  closely  to  those  obtained 
from  a  Klinkenberg  extrapolation. 

Use  of  the  modified  visco-inertial  flow  plot  yielded 
permeability  values  which  were  very  consistent  with  the  Klinken¬ 
berg  extrapolation  results.  The  values  of  the  inertial  re¬ 
sistance  coefficient  obtained  by  the  modified  visco-inertial 
flow  plot  were  found  to  differ  significantly  from  the  standard 
technique  values.  In  view  of  the  fact  that  the  permeability 
results  were  found  to  be  more  consistent,  the  inertial  re¬ 
sistance  coefficient  values  could  be  assumed  to  be  more  re¬ 
liable  . 

The  results  obtained  by  quasilinearization  of  the 
modified  equation  for  small  pressure  gradients  were  found  to 
agree  very  well  with  the  graphical  results.  In  the  present 
work  results  obtained  by  quasilinearization  of  the  equation 
for  high  pressure  gradients  have  been  found  to  compare  very 
well  with  results  obtained  from  the  solution  of  the  equation 
for  low  pressure  gradients.  The  discrepancies  observed  between 
the  experimental  and  predicted  results,  which  may  be  seen  in 
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Figures  (7)  and  (16) ,  have  been  found  to  be  not  very  great  but 
very  similar  for  both  cores. 

Sample  1  was  found  to  have  a  permeability  which 
was  quite  uniform  with  length  and  no  significant  mechanical 
end  effects  were  observed.  Though  core  2  was  found  to  be 
highly  nonuniform,  no  end  effects  of  any  consequence  were  en¬ 
countered  as  there  was  no  noticeable  change  in  the  pressure 
profile  at  inlet  or  outlet  with  a  change  in  the  direction  of 
flow.  Hence,  the  discrepancies  noted  must  be  attributed  to 
deficiencies  of  the  basic  theory  in  describing  the  physical 
model  and  poor  experimental  techniques.  On  the  basis  of  this 
the  merits  of  each  method  may  be  described. 

Adequate  criteria  have  not  been  developed  for  the 
delineation  of  the  regions  of  linearity  on  the  standard 
visco-inertial  plot  and  the  backpressure  plots.  Close  agree¬ 
ment  in  permeability  results  suggests  that  the  modified 
Klinkenberg  and  visco-inertial  flow  plot  technique  is  superior 
to  its  predecessors.  The  shortcomings  of  the  graphical  tech¬ 
nique  is  that  it  is  time  consuming  and  data  must  be  split,  re¬ 
quiring  a  large  number  of  observations  for  good  results. 

Results  obtained  by  the  numerical  solution  of  this 
modified  equation  are  consistent  with  those  obtained  graphical 
ly.  The  modified  form  derived  under  the  condition  of  a  high 
pressure  gradient  is  not  amenable  to  graphical  treatment.  Com 
parison  of  the  results  obtained  using  the  numerical  means 
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with  those  obtained  for  the  modified  form  for  a  small  pres¬ 
sure  gradient  indicates  the  assumption  in  the  derivation 
that  the  mean  pressure  is  constant  is  not  a  serious  one  for 
this  work. 

The  numerical  technique  is  likely  superior  since 
all  data  are  utilized  simultaneously,  and  it  is  simple,  rapid 
and  inexpensive.  Times  of  compilation  and  execution  on  the 
IBM  7040  for  a  typical  set  of  data  were  30  seconds  and  5 
seconds  respectively,  little  more  than  the  time  required 
for  a  linear  least  squares  fit.  The  only  major  drawback  is 
that  convergence  may  not  occur  if  the  initial  guesses  are 
poor.  However,  it  is  thought  that  the  correlations  presented 
for  the  absolute  permeability  and  the  Klinkenberg  b-f actor 
and  the  correlation  developed  for  the  inertial  resistance 
coefficient  provide  adequate  initial  guesses.  The  correla¬ 
tion  proposed  for  the  inertial  resistance  coefficient  has 
a  sounder  basis  than  any  previous  correlation.  Nevertheless, 
it  should  be  used  with  caution  since  it  has  been  shown  that 
the  experimental  values  of  the  three  parameters  which  have 
been  correlated,  the  inertial  resistance  coefficient,  the 
absolute  permeability,  and  the  effective  porosity,  in  many 
cases  are  in  serious  error. 

A  comparison  of  the  predicted  and  observed  results 
obtained  for  gas  flow  in  the  presence  of  an  immobile  liquid 
phase  indicates  that  all  differences  are  very  small.  Hence, 
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the  modified  visco-inertial  flow  equation  is  assumed  to 
adequately  describe  the  phenomenon  observed. 

As  expected,  according  to  viscous  flow  investiga¬ 
tions,  the  effective  permeability  was  found  to  decrease  uni¬ 
formly  with  increasing  liquid  saturation  while  the  Klinkenberg 
b-factor  remained  relatively  constant.  The  inertial  resis¬ 
tance  coefficient  was  found  to  be  a  continuous  function  of 
saturation  which  increased  at  an  increasing  rate  as  the  satura¬ 
tion  was  raised. 


. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


The  conclusions  drawn  as  a  result  of  this  study 

are : 

1)  The  standard  visco-inertial  graphical  technique  is 
adequate  if  there  is  certainty  that  the  data  are  not 
affected  by  molecular  streaming. 

2)  Since  there  is  no  adequate  criterion  for  delineation 
of  the  region  of  linearity  for  the  standard  visco- 
inertial  plot,  where  any  discrepancy  is  anticipated 
or  noted  it  is  proposed  that  the  modified  graphical 
technique  be  used. 

3)  When  access  to  a  digital  computer  is  available,  use  of 
the  numerical  technique  is  superior.  The  method  is 
fairly  simple,  rapid  and  inexpensive,  requiring  no 
data  splitting,  hence  utilizing  all  data  simultaneously. 
Convergence,  however,  may  not  occur  if  poor  initial 
guesses  are  made. 

4)  For  this  work  the  assumption  that  the  mean  pressure  p 
is  not  a  function  of  p  in  deriving  the  modified  form 
of  the  visco-inertial  flow  equation  is  not  a  serious 
one . 

5)  It  has  been  indicated  that  the  data  used  in  evaluating 
the  correlation  of  the  inertial  resistance  coefficient 
with  the  absolute  permeability  and  the  effective  porosity 
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is  doubtful.  Hence,  the  proposed  correlations  should  be 
used  exercising  some  measure  of  caution. 

6)  The  revised  correlation  of  the  inertial  resistance  co¬ 
efficient  with  the  absolute  permeability  and  the  effec¬ 
tive  porosity  predicts  results  better  than  any  previous 
proposal . 

7)  The  description  of  flow  through  a  non-uniform  sample 
lends  itself  to  the  identical  analysis  which  has  been 
applied  in  the  case  of  uniform  samples.  A  single  set 
of  parameters  can  be  used  to  adequately  describe  the 
flow  behavior.  For  all  cases  examined  mechanical  end 
effects  were  found  to  be  negligible. 

8)  The  flow  theory  which  has  been  used  to  describe  single 
phase  gas  flow  is  adequate  for  describing  gas  flow  in 
the  presence  of  an  immobile  liquid  phase. 

9)  The  effective  permeability,  Klinkenberg  b-f actor  and 
inertial  resistance  coefficient  are  continuous  functions 
of  the  saturation  of  an  immobile  liquid,  with  Ffi  in¬ 
creasing  at  an  increasing  rate  with  the  liquid  saturation. 

10)  Serious  errors  may  result  in  using  the  present  dry  core 
correlations  to  predict  the  inertial  resistance  co¬ 
efficient  for  gas  flow  in  the  presence  of  an  immobile 
liquid  saturation. 
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Recommendations 


The  following  recommendations  are  made: 

1)  A  study  of  pressure  profiles  should  be  performed  and  the 
consequences  of  nonuniformities  and  mechanical  end  ef¬ 
fects  of  both  machined  and  fractured  cores  should  be 
carefully  investigated. 

2)  The  applicability  of  slip  theory  should  be  verified  for 
purely  viscous  flow  at  very  high  pressure  conditions. 

3)  Data  reported  by  previous  investigators  should  be  re¬ 
evaluated  by  the  presently  proposed  techniques  to 
account  for  molecular  streaming.  Revisions  to  the  cor¬ 
relations  based  on  this  data  should  be  made  accordingly. 

4)  A  high  quality  set  of  standard  flow  calibrating  devices 
should  be  made  available  to  simplify  the  calibration  of 
all  flow  measuring  devices  over  their  entire  range. 

5)  Investigation  should  be  made  of  the  ramifications  of 
non-isothermal  flow. 

6)  Caution  should  be  exercised  in  using  the  dry  core  cor¬ 
relations  to  predict  the  inertial  resistance  coefficient 
for  gas  flow  in  the  presence  of  a  second  phase.  Using 
the  presently  proposed  techniques,  flow  in  the  presence 
of  an  immobile  liquid  saturation  should  be  studied 
further  so  that  the  relationship  between  the  inertial 
resistance  coefficient  and  the  liquid  saturation  may  be 
established.  Simultaneous  flow  of  liquid  and  gas  under 
visco-inertial  flow  conditions  should  also  be  examined. 
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NOMENCLATURE 

_  2 

area,  L 

2 

Klmkenberg  b-factor,  m/Lt 
collision  diameter,  L 
formation  resistivity  factor 

-2 

viscous  resistance  coefficient,  L 

inertial  resistance  coefficient,  L  ^ 

2 

absolute  permeability,  L 

.  2 

apparent  permeability,  L 

.  2 
effective  permeability  to  gas,  L 

length,  L 

molecular  weight,  m 
Avogadro ' s  Number 
Reynolds  Number 

-3 

concentration  of  molecules  per  unit  volume,  L 
pressure,  m/Lt2 

3 

volumetric  flow  rate,  L  /t 

volume  rate  of  flow  through  a  unit  cross-sectional 
area  of  the  solid  plus  fluid,  L/t 

2  2 

universal  gas  constant,  mL  /t  T 

effective  path  radius,  L 

saturation 

temperature,  T 

velocity,  L/t 

mass  flow  rate,  m/t 

distance  in  flow  direction,  L 


• 

J  ,s*»p  oi  idU»#fli3#q  ovxiosiie 


j  ^Vrf^XJ  luos*  Ofll 
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z 

compressibility  factor 

X 

mean  free  path  length, 

y 

viscosity,  m/Lt 

p 

3 

density,  m/L 

4) 

porosity 

Subscripts 

c 

critical 

r 

reduced 

o 

reference  conditions 

1 

inlet 

2 

outlet 

Superscript 


average  value 
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APPENDIX  A 


A  MODIFIED  VERSION  OF  THE  VISCO-INERTIAL  FLOW  THEORY 


ACCOUNTING  FOR  SLIPPAGE  EFFECTS  AT  HIGH  FLOW  RATES 


A- 2 


A  modified  Forchheimer  equation  has  been  derived 
and,  by  Equation  (24) ,  is 


dp 

dx 


uq 


k  ( 1  +  — ) 
P 


+  FBpq' 


(A-l) 


Multiplying  through  by  p/y  which  is  equal  to  PM/yZRT  by  the 
modified  gas  law,  and  setting  the  constant  mass  flow  rate/area 
pq  equal  to  w/A  results  in 
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Dividing  by  w/A 
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A- 3 


the  differential  equation  may  be  written  as 
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ab  +  py  M 

For  a  constant  viscosity  and  compressibility  evaluated  at  the 

arithmetic  mean  pressure  and  temperature,  integration  using 

(Al) 

items  91.1  and  92.1  and  interchanging  limits  and  sign 

results  in 
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Substituting  Equation  (A-7) 
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APPENDIX  B 

DEVELOPMENT  OF  A  GENERAL  NUMERICAL  PROCEDURE 


TO  SOLVE  THE  MODIFIED  FORCHHEIMER  EQUATIONS 


B-2 


In  general  the  problem  consists  of  solving  a  re¬ 
lationship  of  the  form 


(B-l ) 


fen (k , b , Fg) 
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which  is  non-linear  in  the  three  parameters  k,  b,  and  F  . 

a 

If  the  equation  were  linear  a  unique  solution  for  k,  b,  and 
Fg  would  be  obtained  if  the  number  of  observed  data  points 
n  =  3.  For  n  >  3  a  least  squares  fit  of  the  data  would  re¬ 
sult  in  unbiased  estimations  of  k,  b  and  FD . 

n> 

The  equation  may  be  linearized  through  application 
of  the  Newton-Raphson  approximation 
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where 
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j  =  the  j  data  point  and  j  =  1,  2,  ...  N 
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where 
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E .  =  random  error  associated  with  the  j  data 
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point  and  i^*1  iteration 
Equation  (B-2 )  may  now  be  rewritten  as 
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In  the  event  that  N  =  3  it  is  assumed  that  no  random  error 
exists  or  E?  =  0  for  all  j  and  i,  reducing  the  problem  to 
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a  set  of  three  linear  equations  in  three  unknowns.  In  this 

case  ,  k^ ,  and  may  be  obtained  by  applying  Cramer's  Rule 

i 

and  generating  a  recurrence  relation.  Initial  approximations 

bQ ,  k^,  and  Fg  are  made  to  initiate  calculation  and  iteration 
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is  continued  until  i  =  n  where 
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for  z,  the  convergence  criterion,  chosen  arbitrarily  small. 

If  N  >  3  a  least  squares  criterion  is  applied  and 
the  sum  of  the  squared  errors  for  iteration  i 
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is  minimized  with  respect  to  b^ ,  k. ,  and  Fg  by  setting 
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reduce  to  a  series  of  three  equations,  linear  in  three  unknowns 


for  each  iteration 
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This  system  of  equation  may  now  be  solved  for  each 


iteration  using  a  method  such  as  Cramer's  Rule. 


Initial  approximations  bQ ,  kQ ,  and  Fg  are  made  and 

o 

iteration  proceeds  until  i  =  n  where 


b  -  b 
n  n-1 

k  -  k  , 
n  n-1 

and 

F  -  F 

B  B  , 

n  n-1 

f 

b 

k 

Fn 

n 

n 

B 

n 

for  e,  the  convergence  criterion,  chosen  arbitrarily  small. 
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APPLICATION  OF  THE  NUMERICAL  TECHNIQUE  TO  TWO 


MODIFIED  FORMS  OF  THE  FORCHHEIMER  EQUATION 


C-2 


In  order  to  apply  the  numerical  technique  to  the 
modified  forms  of  the  Forchheimer  Equation  it  is  necessary  in 
each  case  to  evaluate  the  partial  derivatives  with  respect 
to  the  three  constants. 

The  first  modified  equation  is 


f  = 


M(p. 


2. 

p2  ) 


w 


2  y  Z  RTL  ^ 


k(1  +  p> 


+  F 


B 


Ay 


(C-l) 


Now 


■J  _ 


M(P;L2  ~  P22) 
27ZRTL  ^ 


(C-2) 
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1  w 

-  +  F  - 

,  J.  bx  B  - 

k(l  +  =)  Ay 


J 

i-1 


(C-3 ) 
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:i_1(l  +  t>i_1  {=)  3  )  i-1  Ay 


9  f  . 
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-  (=) 
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1  2 

k  ( 1  +  b  4) 


J 
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-(P)j 


ki-i(1  +  bi-i'|)j)2 
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Vi-i 


-1 


k2(l  +  |) 


P  J 


D 

i-1 


-1 


(k .  ,  )  2  (1  +  b.  (i)  j) 

1  -L  1  i.  p 
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For  the  second  modified  version  of  the  Forchheimer 


equation 


f  = 


m(P;L2  "  P22) 
2y~ZRTL  ^ 


Y 


2  2ab  (ab  +  p-^y)  by 

y  -  - 2 - -  ln  -  +  — 

k(pn  “  Po  )  (ab  +  PoY)  kP 


Now  we  have 


Y 

L 

Pi  =  P  + 


(Pl  ’  P2> 


( C—  8 ) 


P2  =  P  - 


<Pl  -  P2) 


and 


(Pf2  -  P22)  =  2P  (P2  “  P2) 


Setting 

and 


NUM.  =  ab  +  p^y 
DEN  =  ab  +  P2y 
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w 
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9  ab  NUM  by 

y - - -  ln  -  +  — 

kP(P1“P2)  DEN  kpj  i_1 
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af 


ab 


-  2ab  NUM 

— -  In  ( - ) 

kp  (p-^-p2 )  den 


ab  a 

— - -  ( - 

kP(P1“P2)  NUM 


a 
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DEN 


j 
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af 
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ak 


j 

i-l 


kp (P1-P2 ) 


( 


"Pi 

k2NUM 


p~  NUM  ab2 

+  — 2 - )  +  In  ( - )  2_ - 

k  DEN  DEN  k  p(p1-p2) 


j 


i-l 


aF 
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9Fb 


j 
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if  2  w 
L  ( 3 — ) 
L  L  Ay 


w  ab 

Y  (2y  —  ~  -3 - 

Ay  kp(p1-p2) 


w 
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NUM 


w 
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The  nomenclature  used  in  the  computer  program  is 


B 

Klinkenberg  b-f actor,  b 

BETA 

inertial  resistance  coefficient,  F 

SD 

D  =  DELPM 

M(PX2  "  P22) 

2UZRTL  w/A 

DCALC 

fnj 

DE  =  DELP 

Pi  -  p2 

K 

absolute  permeability,  k 

p(i) 

9  f  /  9  k 

PM  =  PMAR 

1/P 

Q  ( I) 

9  f  /  9b 

R  <  I ) 

3f/9FB 

V(I) 

W  =  WAV 

w/Ay 

Additional 

nomenclature  used  in  the  computer  output  is 

KA 

apparent  permeability,  k 

a 

PlA 

Pi 

P2A 

P2 

QO 

Qo 

TMK 

T 

C  DATA  CORKELAT ION  USINu  wUA S I L I Nt AR I Z AT  I  ON 
C  SCLUHQN  OF  THE  EQUATION  FOR  HIGH  PRESSURE  GRADIENTS 
C  NSETS  =  NUMBER  OF  DATA  SETS 

C  FMT(J)  =  APPROPRIATE  COMMENT  WHICH  DESCRIBES  DATA  SET 

C_ _ N  =  NUMBER  OF  OBSERVATIONS  IN  THE  DATA  SET _ 

C  K  =  PERMEABILITY  ,  INITIAL  GUESS 

C  b  =  KLiNKENBERG  b-FACTCR  »  INITIAL  GUESS 

C  BETA  =  INERTIAL  RESISTANCE  COEFFICIENT  ,  INITIAL  GUESS 

C  EPSLGN  =  CONVERGENCE  CRITERION 

C  D-DELPM  ,  PM=PMAR,  W=WAV,  DE=OELP 

ICO  FORMAT (IX, 13) 

1C1  FORMAT ( 20 A4 ) 

102  FORMAT (2X,6E12«5) 

1 C 3  FORMAT (2X»4E14*7) 

I C A  FURMAT ( 1 H  0 , 1 9  X  » 17H INITIAL  GUESS  ,  3  (  2  X  ,  E 1 1 . 4  )  ) 

_ 1C6  FQRMAT(20X, 13HI TERAT I  UN  NO. , I  4,3 (2X, Ell . 4 )  ) _ _ 

1 C  7  FORMAT ( 1 H  , 1 8X , 5 ( E 10 . 4  ,  IX  )  , E8 . 2 ) 

1 C 8  FORMAT ( 19X, 36H  PMAR  WAV  DELP  D, 

129HELPM  CCALC  ERROR  ) 

109  FORMAT! 19X.37H  (ATM-1)  (GM/SEC  CM2  CP)  (ATM)  , 

129H  (ATM  SEC/CM2  CP)  ) 

11C  FORMA! (  1H1  ) _ 

113  FORMAT (1H  ) 

114  FORMAT l 1 H  -  ) 

115  FORMAT (  1H-,44X, 1 HK , 1 2 X , 1H 3  ,  1 OX , 4H BE T A ) 

Lib  FORMAT (4lx, 3 bH (DARCIES)  (ATM)  (ATM  SEC2/GM ) ) 

REAL  FMT(20) ,D( 100) ,  PM( 100) ,W( 100) ,DE( 100) , V( 100) ,NUM 
_ REAL  Z(  100)  ,  P ( 100)  ,  U_(  100)  ,R  (100)  ,  T (  100)  , ERROR ( 100) ,K,L _ 

REAL  QCALC( 100) ,X( 100) ,K1 
READ(5f  10ONSETS 
CO  i  JJ  =  1, NSETS 
R E  AC  (  5  ,  10  1  )  (  FMTC  J  )  ,  J  =  1 ,20 

WRITE(6, 110) 

WR1TE(6, 114) 

READ(5,1Q0)N 

C  SETTING  THE  INITIAL  GUESSES 

REAC(5,  102)K,B, BETA, EPSLGN 
DO  2  I  =  1 »  N 

RE AD ( 5 ,  103)U(  I) ,PM(I  )  ,W(  I)  ,DE(  I) 

2  COMT  INUE _  _ 

a  R  I  T  E ( 6 , 101) (FMTI J)  ,  J=l,20) 

WRITE(6, 115) 

WRITE(6,116) 
a  R 1 T  E ( b , 104)K,B,bLTA 
WRI TE (b , 1 13 ) 

_ M  =  0 _ _ _  _ _  _ 

8  S  I GP SQ=0 • 0 
S 1 GPC  =  0  •  0 
S 1 GPR=C  •  0 
S i GGSG=0 • 0 
S IGCR=0.0 
S I GRSU =0 • 0 
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S IGPD^O.O 
SIGGD=0.0 
SI GRu=O.G 
SfGPT=0.0 
S I  Gw  T  =  0 . 0 

_ SIGRT^Q.O _ 

C  SETT  i  NG  THE  ELEMENTS  CF  THE  MATRIX  AND  THE  KNOWN  VECTOR 
DO  3  I- 1  ♦  N 
E-l./K+bLTA*W (  I  j 

NUM=BETA*B*W(I>+(1./PM( I )+DE(I)/2.l*E 
DEN=BETA*B*W(I  )  +  (  l./PM(n-DE<  I  )/2.  )*E 

_ L=k±*2+ B*PM ( I ) *E /K-B**2*6ET A*W ( 1 ) *PM( I) / ( K*OE ( I ) ) *ALOG 

1( ABS( NUM/DEN)  ) 

V(  I  )=E**3/L 
X  (  I  )  =  1 . 0 
Z  (  I  )  =  1 .  u 

P  (  I  )  = (-L*3.*E**2/K**2-E**3*{-2.*E/K**2-2.*B*PM(  I  ) /K**3 
1-B*P M ( 1 ) *  BET A*w ( I  )/K**2-B**2*BETA*W(  I  )*PM( I  ?/ (K*DE(  I ) ) 
1*(  (  -  1 .  /  PM  (  l)-DE(  I  )  /  2  •  )/K**2/NUM+(l./PMU  )  —DEI  I  )  /  2  .  )/K* 
l*2/DEN)+(  A  LOG  (ABS(NUM/DEN)  )  )*b**2*BET  A*W  (  I  )*PM(  I  )/K**2 
1/DE (  I  )  )  )/L**2 

Q (I)  =  <-E**3* ( PM (  I  )  /K*E-B**2*BETA*PM(  I  )  *W (I) / K/DE (1) * < B 
1ETA*W ( I  I / NUM—BE TA* W ( I ) /DEN ) - ( ALOG < ABS ( NUM/OEN ) )  ) *2.*B* 

_ 1BETA*PM(  I)*W(  I  )/K/DE(  I)  )  )/L**2 _ 

RU  )  =  (  L*3.*E**2*W  (  I  )-£**3*(  2.*E*W(  I  )+B*PM(  I  )*W(  I  )  /K-B* 
1*2*BETA*W< I )*PM< I )/K/OE ( I )*< ( 8*W ( I)+W(I )*(1./PM(I )+DE( 
II  i/2.  )  )  /  NUM—  (  B*W(  DHl./Phil  )-DE  (  I  )  /  2  •  )  *  W  {  I)  )  /  D  E  N  )  -  A  L  0 
1G ( ABS ( NUM/DEN ) } * *2 * W ( I ) *P M ( I ) / K/DE ( I) ) )/L**2 
T(I )=V< I ) *X4 I )*Z(  l  >-K*P( I )-B*Q ( I ) -BETA*R( I ) 

_ S_I  GPSw  =  SI  GPSG+P  (  I  )**2 _ 

S  IGPU=S  I GPQ+P  U  >*QU  ) 

SIGPR=SiGPR+P( I ) ^R (  I  ) 

SIGGSQ=SIGGSU+Q< I )**2 
S IGGR-S I GQR  +  U ( I ) *R ( I ) 

SIGkSQ-SIGRSU<-R(  I  )**2 

_ SIGPC=S  IGPQ+P  (  I  )*P<I) _  _ _ _ 

S  I  GwD=S I GOD  +  w i I  ) *D( 1) 

S IGRC=S I GRO  +  R ( I )*D( U 
SICPT^SIGPU-P  (  I  )*T (i) 

SIGQT^SIGQT+Q (I )*T( I) 

sigrt=sigrt*k  (ij^t  (  i  ) 

3  C ONTINUE _ _ _ 

Bl  =  SlGPD—  SIGP T 
B2=SIGQU— S1GQT 
B3=SIGRD-SIGRT 

C  SOLUTION  USING  CRAMER* S  i<ULE 

DDE  T  =  SIGP  SQ*(  SI  GUSU^SIGRSG-SIGQR^Sl  GblR  )-SIGPQ=MSIGPU*S 
1  IGKSw-SIGOR^SlGPR  )  +S  IGPR*  (  S  IGPQ*S  I  GUR-S  I  GQSQ*S  IGPR  ) 
PNDfcT=SiGPQ*(SIGQR*B3-SlGRSQ*62)-SIGQSQ*<SIGPR*83-SIGR 
lSU*bl)  ♦SIGUR*!  S I GP R*B 2-S I GQR*61) 

SNUET=SIGPSG*(S IGgR*B3-S IGRSQ*B2  ) -S  IGPQ*  (  SI  GPR*B3-S  I  GR 
lSQ*Bl)*SIGPR*tSIGPR*B2-$IGQR*Bn 
TNDET^SIGPSQ*( S I GQSQ*B3-S I GQR*B2 ) -S IGPQ* < S IGPQ*B3-S I GQ 
1R*B1 )  +  SIGPRM S I GP Q *62- S I GQSQ*61 )  
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K  1  =  1-  NOE  T/DUE  T 
B  1  =  -  SINCE  T /DDE  T 
BETA  1  =  TNCE  I /DDE  T 

C  A  CHECK  GN  THE  CONVERGENCE  OF  THE  SOLUTION 

i_F (  ABS  (  (  Kl-K)  /K  )  .  L_c  .EPS LON  1  GO  TO  j4  _ 

GO  TG  7 

A  IE  ( ABSU Bl-BI/B) .LE.EPSLCN  )  GO  TO  5 
GO  TG  7 

5  IF  (AGS! ( tfETAl— BETA) /BETA > .LE.EPSLON  )  GO  TO  6 

7  IMABStKl-K)  .GT  .<2.*A8S(K)  )  )  K  1  =  K  +  2  .  *K*K  1  /  ABS  (  K1  ) 

IF ( ABS( Bi -B  )  .GT  .(2  .  *A3S( BJ  )  )  Bl=B+2. *b*B  1/ ABS( B1  ) 

IF! ABS ( BET A  1-BETA) «GT.!2.*ABS( BETA) ) )  BETA  1=BET A  +  2.*B ETA 
1*6ETA1/ABS(B£TA1) 

WRITE (6f 106 ) H, K1 tBl , BETA 1 
K=K  1 
B=B  1 

BETA=BETA1  _  _  __  _ _ 

IF(M.fcg.lOO)  GO  TO  6 
GO  TG  8 

6  WkITE(6,lC6)N*Kl,bl#BETAl 
WRITE(6t 1 14) 

UR  I T  E ( 6  t  1081 

_ WRITE! 6 t 109) _ 

WR1TE!6, 113) 

CO  9  1  =  1  ,N 

IF ( ( ( I+M  ) . EQ.43 ) .OR. { (I+M)«EQ«92) )  GO  TO  11 
GC  TG  10 
il  wR I T E ( 6 , 1 10  ) 

WRITE!  6,114)  _  _  _  _ 

'WR 1 Tb (6 ,101)  (FMT(J)t J  = 1,20) 

WR I TE ( 6 , 1 14 ) 

WRITE (6, 108) 

WRITE! 6, 109) 

WR I TE ( 6 , 1 1 3 ) 

10  OCALC  (  I )  =  V  !  n*x  (  I  )  *  Z  l  I  ) 

ERROR! I  )  =  ABS ! GC A LC! I ) -D ( I  )  ) 

WRITE ! 6, 107) PM!  I ) ,W( I ) ,DE< I  )  ,0! I ) , OCALC (  I ) , ERROR! I ) 

9  CONTINUE 
1  CONTINUE 
STOP 
END 
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APPENDIX  D 


TABULATED  DATA 


0-  2 

FABLE  10 

EXPERIMENTAL 

GAS  FLOW 

DATA 

FOR  CORE 

1 

E  N  C  F  .A  L  E 

DATA  FOR 

SANDSTONE 

SAMPLE 

RUN 

P  1  A 

P2  A 

TMK 

00 

DEL  P  2 

(  ATM  ) 

(ATM) 

(K) 

(CC/SEC) 

( ATM2 ) 

i 

30.84 

2.291 

300. 13 

0.1 566E 

03 

0 • 9456  E 

03 

2 

3  3.14 

2 . 534 

301.01 

0.1775E 

03 

0. 1 092  F 

04 

3 

9.90 

0.981 

300. 16 

0 . 192  5t 

02 

0.9700E 

02 

4 

12.30 

1.026 

300.06 

0.2698E 

02 

0. 1504E 

03 

5 

1 J.  o4 

i  .  06  4 

300.14 

. 0.3613E 

02 

0. L905E 

03 

6 

1  b.  bO 

1.119 

300.14 

0 • 446  1 E 

02 

0.2 39  IE 

03 

7 

17.  b6 

1.216 

300.13 

0. 5635E 

02 

0. 3069E 

03 

8 

21.24 

1 . 432 

300.27 

0. 7675E 

02 

0 • 4489E 

03 

9 

24.  b4 

1.086 

299.92 

0. 1040E 

0  3 

0. 5992E 

03 

IG 

2  7 . 4  8 

1.948 

300. 26 

0. 1270E 

03 

0 . 7  5  L  b  E 

03 

11 

2.9 .4  b 

2.150 

. 3.00 . 28 

0 . I446t 

03 

. (U 8628 E 

03  _ . 

12 

2o.27 

1 . 840 

300.41 

0.1173E 

03 

0 • 6868E 

03 

13 

22.  oO 

1.546 

300. 45 

0.9049E 

02 

0.5175E 

03 

14 

19.26 

1.308 

300.52 

0.oo25E 

02 

0. 3693E 

03 

lb 

1.84 

0.922 

300. 80 

0.6523E 

00 

0. 2542E 

01 

16 

1.12 

0. 922 

300. 80 

0. 1136E 

00 

0. 3981E 

00 

17 

3b .  29 

3. 140 

301.08 

0.2292E 

03 

0. 145  7  E 

Q4 

18 

43.72 

3 . 860 

301. 36 

0.269OE 

03 

0. 1896E 

04 

19 

4  8.6c 

4.628 

301.71 

0.3498E 

03 

0.2342E 

04 

20 

9  2 . 0  Li 

5.196 

301. 84 

0.3  94  6 E 

03 

0.2683E 

04 

21 

6b  .  39 

5.  709 

301. 64 

0.4358E 

03 

0 . 3036E 

04 

22 

64.10 

7.291 

302.23 

0. 5640E 

03 

0 . 405bE 

04 

23 

b9  •  92 

6.463 

302. u3 

0.498  7 E 

03 

0. 3548E 

04 

24 

49.91 

4.  833 

301.69 

0 • 3660E 

03 

0.2468E 

04 

2b 

46.14 

4.263 

3  0 1 . 2  7 

0.3194E 

03 

0.2110E 

04 

2o 

4  0  .  8  0 

3.496 

300.91 

0 . 2  5  7  6  E 

03 

0. 1652E 

04 

3b.4b 

3.200 

300.56 

0.2001E 

03 

0. 1247E 

04 

28 

3b.  49 

3.205 

300.40 

0. 1996E 

03 

0. 1249E 

04 

49 

3  3.07 

2.540 

300.24 

0. 1770E 

0  3 

0. 1087E 

04 

3C 

21.16 

1.427 

300.09 

0. 7  8  8  9  E 

0  2 

0.4459E 

03 

31 

b  3 . 6  6 

5 . 4  b  o 

299.76 

0.4  L62E 

03 

0. 2850E 

04 

32 

71.59 

8 . 73  3 

300.25 

0 . 6  8  4  9  E 

03 

0. 5U49E 

04 

33 

od.90 

6 . 206 

300. 01 

0.6434E 

03 

0.4680E 

04 

34 

66 *  74 

7.814 

299.90 

0.6 100E 

03 

0. 4393E 

04 

3b 

04.26 

7.  355 

294. 72 

0.5721b 

03 

0.4075b 

04 

36 

O  1  •  86 

6  •  904 

29v.o4 

0.53 5 6 E 

03 

0. 3779 E 

04 

37 

57.71 

6.173 

299. 51 

0.4759E 

03 

0. 3292E 

04 

38 

b4 .  b  1 

5.620 

299.47 

0.4324t 

03 

0^29 4J E 

314 

39 

8  •  b  0 

0 . 960 

300. 20 

0. 143 3E 

02 

0. 7138E 

02 

4  C 

1.31 

0.920 

300. 53 

0 • 2  366E 

00 

0.8728E 

00 

41 

1.63 

0.920 

300.54 

0.4748b 

00 

0.  1822E 

01 

42 

2.01 

0.921 

300.6b 

0. 7976E 

00 

0. 3 1 9 1 E 

01 

43 

2.7b 

0.923 

300. 56 

0. 1 58bt 

01 

0.6718E 

01 

44 

4.33 

0.9  26 

30u.o2 

0. 3893E 

01 

0. 1 788E 

_ 
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TABLE  10 

EXPERIMENTAL  GAS  FLOW  OATA  FOR  CORE  1 

t NCR AC  E 

DAT  A  FOR  SANDSTONE  SAMPLE 

RUN 

PiA 

P2  A 

TMK 

00 

DELP2 

(  ATM  ) 

(ATM  ) 

(  K  1 

(CC/SEC ) 

( ATM2 ) 

43 

S  •  6  3 

0.934 

300. 34 

0  •  2  7  d  S  E 

01 

0. 1243E 

02 

46 

i  .  3u 

0.923 

300.6b 

0.3729E 

00 

0. 1 402  E 

01 

47 

1.40 

0.922 

300. 71 

0.3024E 

00 

0. 1 122E 

01 

4b 

2.42 

0.926 

300. 36 

0. 1215E 

01 

0.3023E 

01 
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10 
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D-  4 

T  AbL  c  ID 

EXPERIMENTAL  GA5  FLOW  DATA  FOR  CORE  I 

ENCFACE  DA 1  A  FUR  SANDSTONE  SAMPLE 


RUN 

DELP 

PMAR 

KA 

DELPM  WAV 

(  ATM  ) 

(ATM-1 ) 

( MD) ( ATM 

SEC /C  M2  CP) ( DM/ S  EC  CM2  CP) 

i 

0 . 2  8  54i_ 

02 

0.O037E-01 

0.2161E 

01 

0. 4o23t 

03 

0. 1630E 

00 

2 

D . 5  J6lc 

02 

0*30  06  E  —  0 1 

0.2133E 

01 

0.4688E 

03 

0. L840E 

00 

3 

0 . 89  17E 

01 

0.1B38E  00 

0  •  2  b  6  7  E 

01 

0.3896E 

03 

J.2023E- 

■01 

4 

0. 1128E 

02 

O.lbOOE  CO 

0.2494E 

01 

0.4010E 

03 

0.3043E- 

■01 

b 

C . 12  /8E 

u2 

0.1342E  00 

0 . 2  4  b  6  E 

01 

0.40  7c  E 

03 

0.3790E- 

■01 

6 

0 . i438E 

02 

0.1203E  00 

u. 2418E 

01 

0 . 4 1  3  6  E 

03 

0.46  76E- 

•01 

7 

0 . 16  34L 

02 

0 . 1 uo be  00 

0.2381E 

01 

0.41 99E 

03 

0.5902E- 

Qi 

8 

G. 19  60b 

02 

0. 8d23£~0i 

Q.ZZ6QL 

01 

0. 438bt 

u3 

0.8230F- 

■01 

9 

0. 228bE 

02 

0. /o27E-01 

0.2234E 

01 

0. 4437E 

03 

0. 1086 E 

uO 

10 

0 . 2  j  b  jL 

g2 

0. 6/96E-0I 

0.2203E 

01 

■  J  •  4  5  3  8 1 

03 

0. 1323E 

00 

U 

0  •  ^  7  3  C  t 

02 

0.6329E-01 

0.218oE 

ol 

0.4574E 

03 

0. 1505E 

00 

12 

0 . 2443E 

02 

0. 7114E-01 

0.2227E 

01 

J.4490E 

03 

0. 1223E 

00 

1J 

0.2i2bE 

02 

0.8215E-C1 

0.227 IE 

01 

0. 4392E 

03 

0. 9443b- 

•01 

14 

0 . 179  bb 

02 

0 . 9  72  3E  —  U 1 

0. 2333  E 

01 

0.4286E 

03 

0 • o9 26E- 

01 

lb 

C  .  9  1  9  9  E 

00 

0.  /2  3  7b  00 

0 . 332  IE 

01 

0. 30  1 1  E 

03 

0.6869E- 

■03 

16 

0  •  1 9  b  2  L 

00 

0.9808E  UO 

0. 3691 E 

01 

J.2709E 

03 

0. 1 196E- 

•03 

17 

J  .  3  ‘J  l  JL 

02 

0 . 482  7 E—  6 1 

U.2071E 

01 

0 •  482  7 1 

03 

0.23  I0E 

00 

18 

0 . 3  9  8  4  E 

02 

0. 4202E-C1 

0 . 2  0 1  8 1 

01 

0.4934E 

03 

0.2983E 

00 

19 

0.4399E 

02 

0. 3 7  bo E~ 01 

0. 19d3E 

01 

0. 3042E 

03 

0.3391E 

00 

20 

0  .  Hu  bo t 

02 

0. 3493E-C1 

0. 1 9  b9 L 

01 

J. blOob 

03 

u • 4041 E 

00 

21 

0 • 4  968E 

02 

0*32  73E-01 

0. 1 9 1 3  E 

01 

0. 3228E 

03 

0.4437E 

00 

22 

O.bodlE 

02 

0 . 2801 E- 01 

O.i 869E 

01 

0  .  b  3  b  0  E 

03 

0.3730E 

00 

23 

0. b  3  4 3  E 

02 

0.30LZE-CI 

0. 1882E 

01 

0.3313c 

03 

O.3082E 

00 

24 

0 . 4  bCdE 

02 

0 • 3ob3  E~  0 1 

0. 1971E 

01 

J.  b0  /4E 

03 

0.3753E 

00 

2b 

C  •  h  1  8  7  E 

02 

0. 3968E-01 

0.2002E 

01 

0 • 49  94  E 

03 

0.3287E 

00 

26 

0.  W30E 

02 

u .43 13l-01 

0. 203JL 

01 

J.4871E 

03 

0.2661E 

00 

27 

0 • 3  2  25E 

02 

0. bi 74E-C1 

O.2103E 

01 

0.4753E 

03 

0.2074E 

00 

28 

C . b  22  dE 

u2 

0 . b lo9c-C l 

0.2  092  E 

01 

0 . 4  7  7  9  l 

03 

0.20  7 OE 

00 

29 

o  •  3  0  5  3  L 

02 

0 • b  o 1 6E~  01 

0  •  2  1  c.  I L 

01 

0.4702E 

68 

0. 1839E 

00 

30 

0. 1974E 

02 

0. 88b3c-01 

0.2297E 

01 

0.4332E 

03 

0.82 *+92- 

01 

31 

0.482  IE 

02 

0. 3384E-01 

0. 1922 E 

01 

0.  32 02 E 

03 

0 . 4  2  8  0  E 

00 

32 

0 • j2 86 t 

02 

0.2490E-C1 

0 . 1  8  1 0  E 

01 

0.  b 323E 

03 

0 . 69  b  9c 

00 

33 

0.6C70E 

02 

0.2594E-C1 

0. 1829E 

01 

0  •  b  4  6  9  E 

03 

0 . 6bb2E 

00 

34 

0 • bd  93E 

02 

0. 2683t-0l 

U • 1 843  E 

01 

0. 3423E 

03 

0.6223E 

00 

3b 

C . 3O90L 

02 

J . 2  793E-01 

u  .  1 8  3  9E 

01 

).  5379E 

03 

0. b  848c 

00 

36 

0  .  b  4  9  6 1 

02 

0.2909E-G1 

0 •  i 8  73 1 

01 

0 . 33  39 1 

03 

0. 5484E 

00 

37 

C . b i b«+E 

02 

0.313  it- 01 

0. 1904E 

01 

0.5232E 

03 

0.4886E 

00 

38 

C.4889E 

02 

0.332oE-0i 

0. 19  33 E 

01 

0. bl 72l 

03 

0 • 4448  E 

JO 

39 

0  .  / b4 3E 

01 

0 *  2  1  1 4 E  CO 

0.2396E 

01 

0. 38b2t 

03 

0. I507E- 

01 

40 

0. 3912E 

00 

0.8964E  00 

0 . 3  b03E 

01 

0.2833b 

03 

0.2494b- 

■03 

41 

6 •  /  1  JCL 

00 

0.7832E  Go 

U  •  3  3  6  7  E 

01 

0.29  70E 

03 

0. 5004E- 

■u  3 

42 

D . 1 J89L 

01 

0. od2bc  DO 

0. 3232E 

01 

0. 3094E 

03 

0. 84 J1E- 

■03 

43 

C . 1 o28E 

01 

D  •  3443 E  00 

0  •  30b 1 E 

01 

0.3277E 

03 

0.1670E- 

•02 

44 
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TABLE  10 

EXPERIMENTAL  GAS  FLOW  DATA  FOR  CORE  1 

cl^Ll  ACL 

DATA  FOR  SANOSTONE  SAMPLE 

ROM  OLLP 

(  ATM) 

P  MAR  KA  OELPM  WAV 

(ATM-1)  (MO)  (AIM  SEC/C  M2  CP)(GM/SEC  CM2  CP) 

A  3 

U  .  2  / 1AL 
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TAdLt  11 

EXPtRI MENT AL  GAS  FLOW  DATA 

FOR  CORE  1 

DATA  FOR 

SANDSTONE  8 AM  PLt  FROM  SECTION  1 

RON  PiA 

(ATM) 

P2  a  tmk  go 

(ATM)  (K)  (CC/SEC) 

DELP2 
( ATM2 ) 

1 

30 .  64 

26.481 

300. 13 

0.1 866E 

03 

0.2496E 

03 

2 

33.14 

26.4/9 

301.01 

0. 1775E 

03 

0. 2872F 

03 

3 

9.90 

8.494 

300. 16 

0. 1925E 

02 

0 . 2  3  8 1 E 

02 
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12.30 

10.670 

30u . 13 

0.2898E 

02 

Q.3969E 
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lA.dk 

11.880 

100.  14 

0 . 3  o  1 3  fc 

02 

Q . 3034E 
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0  •  ‘+  4  o  1  £ 
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0.O253E 
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2  1 . 2  A 
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0  .  /  b  /  3  E 
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0. 1213P 

0  3 
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299. 92 

0 • 1 040E 

03 
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03 
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0.  12  70t 
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0. 1989E 

03 
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2  3..  4  3 . 

23.319 . 

300.28 

0 . 1 44ol 

03 

0.2264E 

03 

12 

2  6  •  2  7 

22. 39  7 

300.41 ' 

0.1173E 

03 

0. 1795E 

03 

13 

22.80 

19. 686 

300.43 

0.9049E 

02 

0. 1363E 

03 
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T  AdLE  11 

EXPERIMENTAL  GAS 

FLOW  DATA  I 

FOR 

CORE  1 

DATA  FOR 

SANDSTONE  SAMPLE  PROM  SECTION  1 

RUN 

OELP 

PMAk 

KA 

DEL  PM 

WAV 

i  ATM) 

( ATM-1 ) 

i  ML)  i  (  ATM 

SEC /C M2  l 

CP)  ( 

GM/SEC  CM2  CP) 

1 

Li  .  4  3  3  3E 

Cl 

0. 3489E-01 

0. 1952E 

01 

0. 3123E 

03 

0. 1611E  00 

2 

0  «  4 6 6  1 1 

01 

0-32  46E— 01 

0-1936E 

Ol 

0 . 3 1 06  L 

03 

0.1817E  00 
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0.14036 

01 

u* 10676  00 

0  •  l  2  8  3  E 

01 

0.4381E 

03 

0.2017E-01 

4 

0. 1  7336 

01 

0.87436-01 

0.2236E 

01 

0 • 4468c 

03 

0.3029E-01 

. .  . 5. . 

C*  1 96  5 6 

ul 

0. 77756-01 

0.2194E 

01 

0.4339c 

03 

0 . 3  7  72E-01 
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Q  .  2 1  o9£ 

01 

0.O936E-01 

0.2191E 

01 

J.4563E 

03 

0 .46d0£-0 1 
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0.246/6 

01 

0*  6124E-01 

0. 2133l 

01 

0*  464 6 E 

03 

0 . 3  8o5  E-0 1 
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0 .30796 

0 1 

0. 50776-01 

U.2008E 

01 
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0. 3917E-01 
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01 
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12 

EXPERIMENTAL  GAb  FLOW  DATA  FOR  COKE  1 

LATA 

FUR 

SANuSTGNE  SAMPLE  FROM  SECTION  2 

RUN 

P  LA 

P2A 

TMK 

yo 

DELP2 

(  ATM) 

(ATM) 

(K) 

(CC/SEC ) 

( ATM2 ) 

1 

26.48 

22.228 

300. 18 

0.1 5ooE 

03 

0 . 2  0  7  l  E 

03 

2 

28.48 

2  4. 1 0  1 

301.01 

0. 17  73E 

03 

0.2299E 

03 

3 

8. 49 

7.103 

300. lb 

0. 1923E 

02 

0. 2170E 

02 

A 

10.37 

8. 831 

300.  13 

0.2648E 

02 

0.3373E 

02 

5 

11.86 

9.94  i 

300. 14 

0 . 3  6 i 3E 

02 

0. 4218E 

02 

6 

13.33 

11 . 138 

300. 14 

0 • 446  IE 

02 

0.5374E 

02 

7 

13. 09 

12.633 

300. 13 

0 . 3b35E 

02 

0.o82  IE 

02 

t> 

id.  16 

18.200 

300.27 

0. 7873E 

02 

0.9862E 

02 

9 

21.  12 

17.677 

294.92 

0. 1040E 

03 

0. 1334E 

03 

10 

2  3.  3  8 

19. 739 

3  00 • 2b 

0. 12  7JE 

03 

0.1658E 

03 

U 

23.32 

21. 229 

300.28 

0 . 144oE 

03 

0. 1904E 

03 

12 

22. 6U 

18.913 

300.41 

0.1 l 7  3E 

03 

0. 1328E 

03 

13 

19.  39 

16.403 

300.43 

0.9049E 

02 

0. 1  143E 

03 
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TABLE  12 

EAPcRI MENTAL  GAO  FLOW  DATA 

FOR 

CORE  1 

CAT A  FUR 

SANuSTCNE  SAMPLE  FROM 

5ECTIQi\  2 

RUR 

UcLP 

PMAR 

KA 

DELPM 

WAV 

(  ATM) 

( ATM-1) 

(MO) (ATM 

SEC /C M2  i 

CP) (6M/3EC  CM2  CP) 

1 

0 . 4  2  5  3  fc 

01 

0.41 06 1—  0 1 

0.2028E 

01 

0. 4932E 

03 

0 .  16  1  BE  00 

2 

0,43  72E 

01 

o.3  8  03  h—  0 1 

0.2084E 

01 

0.479 8b 

03 

0.1825E  00 

3 

o  •  13  42c 

01 

0.1262c  00 

0.2347E 

01 

0.4261E 

03 

0.20  19E-0 1 

4 

0.17 B9b 

01 

O.iOilE  00 

0.22  77 E 

01 

0.4392E 

03 

0 • 3034E— 0 1 

5 

0 .15336 

01 

0.9 163c- 01 

u  .  2  2  7  1 1 

01 

0.4402c 

03 

0.3779E-01 

6 

0 . 2 1 9oc 

01 

C. 81 73E-01 

0.2203c 

ul 

0 . 453  8E 

03 

0.4660E-01 

7 

C.2460E 

01 

0. 7212c-C1 

0.2196b 

01 

0.4554E 

03 

0.58  78E-01 

a 

0 . 2  9  5  7 b 

01 

0 . 3996C-01 

0.2129E 

01 

0. 4697 E 

03 

0.81 90E-01 

9 

C • 3  44Cb 

01 

0. 515b E— Cl 

0 . 2  0  7  7  E 

01 

0.4813E 

03 

0.1079E  00 

10 

0 . 3  324t 

01 

0 • 40 14c— 01 

0.2052E 

01 

0.48  74E 

03 

0.1315E  00 

11 

_  g.hGsoe 

gi 

0.42 97 £-01 

.  0.2 036c. 

01 

0.4910c 

03 

0. 1494E  00 

12 

0 . 368 it: 

01 

0.4818E-C1 

0.2055E 

01 

0.4865b 

03 

0.1215E  00 

13 

0 . 3  1  8  1 L 

01 

0. 5557E-01 

0.21 12E 

01 

0.4735E 

03 

0.9395E-01 

Lit 

0  .  2  7  Q 5  0 

ill 

..-.Q 6.5  /I£-01 
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EXPERIMENTAL  GAS  FLOW  DATA  FOR  CORE  1 

GAT  A 

FUR 

SANDSTONE  SAMPLE  FROM  SECTION  3 

RUN 

PI  A 

P2A 

T  MR 

go 

0ELP2 

i  ATM) 

(ATM) 

(K) 

(CC/SEC) 

(  A  T  M  2  ) 

1 

22.23 

17.655 

300. 13 

0.1 566E 

03 

0. 1831E 

03 

2 

24.  11 

19.150 

3  0 1 . 0 1 

0. 1775E 

0  3 

0. 2144F 

03 

3 

7.10 

3.585 

300. 16 

0. 1925E 

02 

0. 1925E 

02 

4 

8.83 

6.938 

300.15 

0.2898E 

02 

0 . 2  98oE 

02 

. . . .  5l . 

9.95 

/  .822 . 

300.14 

0...  36  UE 

02  ... 

0. 377o£ 

02 

6 

11.14 

6.783 

300. 14 

0.4461E 

02 

0. 4o90E 

02 

7 

12.63 

9.966 

300.13 

0 . 56  3  3L 

02 

0.6032E 

02 

8 

15.20 

12.024 

300.27 

0.78/5 E 

02 

0. 8647E 

02 

9 

17.68 

13. 969 

299.92 

0. 1U40E 

03 

0. 1 168E 

03 

IC 

19.76 

1  3 . 65o 

300.26 

0.12  JOE 

03 

0. 1453E 

03 

_ . i  JL 

10.330 . 

300.^8 

0. 1446E 

03 

0. 1674E 

03 

12 

18.52 

14.993 

300.41 

0.1173E 

03 

0. 1330E 

03 

13 

16.40 

12.994 

300. 45 

0.9049E 

02 

0. 1003E 

03 

L4 

13.85 

10.944 

3011.2^2 

0.6625E 

02 

Q.  Z2L3£ 
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EXPERIMENTAL  GAS  FLOW  OATA  FOR  CORE  1 

OA  TA 

FbR 

bANGS  TONE  SAMPLE  pRuM  ScCTIGN  3 

RUN 

JELP 

PM  AR 

KA 

DEL  PM 

WAV 

{  ATM) 

( A  TM  -1  ) 

(MO) (AIM 

SEC/ CM2  i 

CP) (GM/SEC  CM2  CP) 

1 

0 • 43  93  E 

01 

0.3017E-01 

u.  2283c 

01 

0.4376c 

03 

0.1523E  00 

2 

0  .  4  9  )  7 1 

01 

C.4g24c-C1 

0.2225E 

01 

0.4493b 

03 

0.1834E  00 

3 

G. 13  1  7E 

01 

0  •  1  3  7  6  E  CO 

0.2643E 

01 

0. 3784b 

03 

0 • 2022E-0 1 

4 

C. 1893E 

01 

0.  12b8E  00 

0 . 2  368 1 

01 

0. 3893E 

03 

0 • 3  039E— 0 1 

5 

0-2 123c 

01 

0.1126E  GO 

0.2534E 

01 

0. 3947E 

03 

0.3  /86E-01 

6 

0*2  334t 

01 

C.1004E  GO 

0.252QE 

01 

0. 39o8b 

03 

0  •  46b9E-0 1 

7 

0. 2669E 

01 

0. 8349E-01 

0 . 2  4  78 1 

01 

0 . 4033t 

03 

0.3892E-01 

a 

0.31  /6E 

01 

J.  /346E-v;i 

0 . 2  42z  F 

01 

0.4128E 

03 

0.82  13E-0 1 

9 

o • 3o8oh 

01 

0.63 16E-0 1 

0.2367E 

01 

0 . 4  2  2  3  E 

03 

0.1083E  00 

10 

u  .  4 i 0  3E 

01 

0 . 3b47  c— Cl 

0.2332E 

01 

0.4287E 

03 

0.1320E  00 

li 

G  •  4  3  9  9  E 

01 

0.3255E-01 

0.2307E 

01 

0  *  4333b 

03 

0.1500E  00 

i  2 

0.3923 E 

01 

G.5898£-Gi 

0.2384E 

01 

0  •  4-248 1 

03 

0.1219E  00 

13 

0. 34  lit 

Gi 

0. od03t-01 

0 . 2  403 1 

01 

0 . 4 1 39 t 

03 

0. 9424 E— 01 

Lit 

v  .  2  9  0  9  £ 

1LL 

o . 3J66L- 0 1 

0. 2444b 

0.4092E 

03 

Q.6913E-01 
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TABLE 

14 

EXPERIMENTAL  BAS  FLOW  DATA  FOR  CURE  1 

OAT  A 

FuR 

SANu STONE  5  AM  PL  L  FROM  SECTION  4 

RUN 

P  1 A 

P2  A 

TMK 

QG 

DELP2 

i  ATM  J 

(ATM) 

IK) 

(CC/oEC) 

( AT M2 ) 

1 

W.b3 

2.291 

3  uO .  13 

0 . 13b6E 

03 

0. 3057E 

03 

2 

14.13 

^  .  334 

301.01 

0. 1773E 

03 

0. 3b03E 

03 

3 

3.34 

0.981 

300.  Id 

0 . 1923E 

02 

0. 302  3 E 

02 

4 

o  •  94 

1 . 026 

3  CO .  13 

0.2898c 

02 

0.4708E 

02 

5 

7.82 

1. 0b4 

300. 14 

0  •  3  6  1 3  E 

02 

0 • bOObE 

02 

6 

3.7  b 

1.119 

300.  14 

0 • 446  it 

02 

0.  7  5 9 0 E 

02 

7 

9.97 

1.216 

300.  13 

0 . 3b  33t 

02 

0.  97  84 1 

02 

8 

12.02 

1.432 

300.27 

0 . 7  8  7  9 1 

02 

0. 1425E 

03 

9 

13.99 

1  .  odb 

299.92 

0. 1040E 

03 

0. 1929E 

03 

10 

1  3.00 

1.948 

3  Ou. 26 

0. 1270E 

03 

0.2413E 

03 

ii 

i  u  »  o  3 

2.130 

300.^8 

0 • 1 44oE 

03 

0 . 2  786  fc 

03 

12 

14.99 

i  .  340 

300.41 

0.1173E 

03 

0. 2214E 

03 

13 

12.94 

1 . 346 

300.43 

0.9049E 

02 

0.1 bb4E 

03 

14 

10.94 

i  .30 8 

300. 32 

0.6623E 

02 

0 . 1 1 8 i£ . 

03 
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rAdLt  14 

EXPERI MENTAL  GAS  FlOW  OATA  FOR  CORF  1 

C  A I A  FOR 

SANuSTUNE  SAMPLE  FROM  SECTION  4 

KUR  CELP 

(  ATM ) 

PMAR  KA  DELPM  WAV 

(AIM-1)  (MG)  (ATM  SEC /C M2  CPHGM/SEC  CM2  CP) 

1 

0. 1 339c 

02 

0. 1  J04c 

00 

0.23  o4E 

01 

0.4230E 

03 

0.1640E  00 

.2 

0 .  1 6o  2  £ 

02 

L.9224E- 

-01 

0-2286F 

01 

0.4373F 

03 

0-1R32F  OO 

3 

0 . 4  o  C4 1 

01 

U • 3u46E 

00 

0.2924E 

01 

0. 3420E 

03 

0.202  7E-0 1 

4 

0.3912E 

01 

0. 23  1  IE 

00 

0.2628E 

01 

0. 3536E 

03 

0.3050E-01 

. 3 

. c . o 

01 

u • 22  3  i  t 

00 . . 

0*2 7o4E . 

0.1 

0.361 8 1 

.03 

0  •  3.80JLE-.0 1 

6 

C  •  7 6  o  4 1 

01 

0  •  2  02  0 1 

00 

G.2702E 

01 

0*  3 70 IE 

03 

0. 46 90 E- 01 

7 

0. 6  7  496 

01 

0.1  789E 

CO 

C.2649E 

01 

0.3775E 

03 

0. 3922E-01 

8 

0. 1U59E 

02 

0 .  1  4  8^E 

00 

0.2343E 

01 

0. 3929E 

03 

0.  b  9  6 4 1 —  0  1 

9 

0. 1230E 

02 

0. 12  76E 

00 

0.2480E 

01 

0.4032E 

03 

0.1091E  00 

10 

U. 13  71E 

02 

0  .  1 1 36 E 

00 

C.2426E 

01 

0  •  4 1 1 8E 

03 

0. i 331E  00 

u 

0  .  1 4qcL 

.02  

0.  10  34  £ 

0  0  „ 

0*  2  395 £ 

01 

0*41  78  E, 

03 

0. 1514E  00 

12 

0. 131:>E 

02 

0. 1188 E 

00 

0.2446E 

01 

0.4089E 

03 

0.1229b  00 

13 

C.  L143E 

02 

0. 13  76E 

00 

0.2308E 

01 

0. 3988E 

03 

0*  94  8  7E-01 

L4 

0 . 9  6  3  j  E 

01 

0 . 1 o32 £ 

Alii.... 

-U*2.3JiZ£- 

01 

Q.3863.E-. 

Q3 . 

.  0 . 6952E-.0] 
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TABLE  15 

EXPERIMENTAL  GAS  E LOW  DATA  FOR  CORE  2 

cNUF ACE 

DATA  FOR  SANDSTONE  SAMPLE 

RUN 

PLA 

P2  A 

TMK 

DO 

l )ELP2 

(  ATM  ) 

l  ATM  ) 

(K) 

(CC/SEC) 

(  A  T  M  2  1 

1 

1.70 

0.949 

300.23 

0 • 8468E 

01 

0. 19826 

01 

2 

1.38 

0 . 946 

300.29 

0.6 84 IE 

01 

0. 1 396F 

01 

3 

1.43 

0.941 

300. 34 

0. 3343E 

01 

0.  1218E 

01 

4 

1.36 

0.939 

300. 31 

0 .42 1 3E 

01 

0.9719E 

00 

3 

i  •  j>3 

0.93  7 

300.37 

0 .4022 E 

01 

0.897  QE 

00 

6 

1.29 

0.93  3 

300. 34 

0.3528E 

01 

0. 7858E 

00 

7 

1.26 

0  •  934 

300. 40 

0.3149E 

01 

0.7153E 

00 

8 

1.21 

0.92  7 

300.23 

0.2713E 

0  1 

0.  o  9  6  4  F 

00 

9 

1.  16 

0.92  6 

300.31 

0.2213E 

01 

0.4843E 

00 

10 

1.11 

0.914 

300.26 

0. 1799E 

01 

0.3963E 

00 

. . IX  _ 

1.08 

0 . 91  3 

.  300.27 

Q. 1357E 

01 

0. 3432E 

00 

12 

1.  7  4 

0.93  7 

299. 98 

0.9208E 

01 

0.2156E 

01 

13 

1.70 

0 . 93o 

299. 98 

0. 8704E 

01 

0.2008E 

01 

14 

1.62 

0.9  33 

3  QO . 0  1 

0.74 9  9 F 

0  1 

0.  1  7  44  F 

01 

13 

1. 30 

0.929 

299.97 

0 . 6034E 

01 

0. 1378E 

01 

16 

1.39 

0.923 

299.90 

0.4782E 

01 

0.  1077E 

01 

17 

. 1.31 

0.923 

299. 89 

0 . 3  8  2  6  E 

01 

0.8  348  E 

00 

18 

1*34 

0.918 

299. 39 

0.4206E 

01 

0.941  IE 

00 

19 

1.28 

0.926 

300. 10 

0.3 3  JOE 

0  1 

0 . 7668E 

00 

20 

1.22 

0.923 

300.00 

o.283Jt: 

01 

0.o269£ 

00 

21 

4.o3 

1.386 

300.26 

0.71 7  IE 

Oz 

0.  1 952  E 

02 

22 

1.  7 3 

0.937 

299.39 

0.92  75E 

01 

0.2180E 

01 

23 

2.17 

0.939 

299. 71 

0. 1362E 

02 

0.3  796  E 

01 

24 

2.39 

0.991 

299.o2 

0. 228JE 

02 

0.5710E 

01 

23 

2.80 

1.010 

299.31 

0.269JE 

02 

0.6  793E 

01 

26 

2.38 

0.973 

299.61 

0. 1909E 

02 

0.4709E 

01 

2  7 

1.96 

0.946 

299.63 

0. 1233E 

02 

0.2944E 

01 

28 

2.98 

1.031 

300.46 

0. 3087E 

02 

0. 783 1 E 

01 

29 

i.Oo 

0.928 

300. 31 

0. 1243E 

0.1 

0.2705E 

00 

30 

i  .  03 

0.928 

300.28 

0. 1070E 

01 

0.231  7E 

00 

31 

1 .  u2 

0. 927 

300. 30 

0. 8746 E 

00 

0. 1885E 

00 

32 

i.Ol 

0.927 

300. 06 

0. 7  8  3oE 

00 

0.  1 1>88  E 

00 

33 

1.00 

0.92  7 

300.20 

0.6991E 

00 

0. 1493E 

00 

34 

1.00 

0. 933 

300. 24 

0 • 30 3  3E 

00 

0. 1210E 

00 

33 

0.98 

0.9_>^ 

2 

0.3992E 

00 

0.  3 0 5 6 E ■ 

-01 

36 

C.  97 

0.931 

300. 12 

0  •  3  0  3  3  E 

00 

0.O598E- 

-01 

37 

9.02 

300 

299.97 

0.2310E 

03 

0. 7051E 

02 

38 

10. 74 

*+.279 

3  00 . 0  3 

0.308  7E 

03 

0.9697E 

02 

39 

11.  73 

4.903 

299.67 

0. 33  77E 

03 

0. 1 1 35E 

03 

40 

13.10 

5.  736 

300.26 

0.4243E 

03 

0. 1384E 

03 

41 

2  8. 72 

17.475 

299. 96 

0. 1331E 

04 

0.3196E 

03 

42 

1  /.  83 

9.022 

300.04 

0 • 6798E 

03 

0.2374 E 

03 

4  3 

2  7.41 

16.413 

300. it 

0. 1243E 

04 

0 . 482  IE 

03 

44 

26.0  / 

13.304 

300.41 

0. 1  lc>2E 

04 

0. 443+E 
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TAoLE  15 

EXPtRIMEN 1 AL 

Gas  flow  I 

DATA 

FOR  CORE 

2 

t N  6  P AG  6 

DATA  PGR 

SANDSTONE 

SAMPLE 

RUN 

P  1A 

P2A 

TMK 

QO 

06LP2 

{  ATM) 

(ATM) 

IK) 

(CC/StC) 

( ATM2 ) 

45 

2  •+  •  8  4 

14.358 

299.30 

0 . 10  86c 

04 

0.41116 

03 

4  6 

2  4.03 

13.  7  i  2 

299. 80 

0.10386 

04 

0.38936 

03 

47 

14.33 

6 . 492 

300. 79 

0 • 462  5E 

03 

0. 1 o  3  2  6 

03 

48 

19.28 

9. 983 

300. 53 

0.75 7  IE 

03 

0.2720E 

03 

49 

2  0.69 

1  1.031 

30u, 33 

0.839/E 

03 

0. 3063E 

0  3 

50 

19.97 

10.507 

300.  19 

0. 7995E 

03 

0.2885E 

03 

51 

18.47 

9.405 

300.06 

U. 7126E 

03 

0 . 2  o  2  6  E 

03 

52 

23.02 

12.909 

299.54 

0.9741E 

03 

0.3635E 

03 

53 

i  8.  79 

6.  139 

300. Go 

0 . 45  66  E 

03 

0. 1524E 

03 

54 

12.09 

5.050 

300.17 

0. 3710E 

03 

0. 1207E 

03 

55 

i  1.07 

4.39  7 

299.93 

0.3193E 

03 

0. 10316 

03 

56 

9.44 

3  •  7o4 

300.  11 

0. 269  IE 

03 

0.8455E 

02 

57 

8.34 

1  •  08  4 

300. 39 

0. 3854E 

02 

0.99896 

01 

58 

3.92 

1.203 

3  CO . 4U 

0.5248E 

02 

0. 1392E 

02 

59 

4.45 

1.339 

300.21 

0 . 66  bOE 

02 

0. 1802E 

02 

60 

4.84 

1.433 

300.07 

0. 779oE 

02 

0.2133E 

02 

61 

6.4  1 

1.693 

300. 15 

0 . 9  5  3  2  E 

02 

0 . 2  o  5  3  E 

02 

62 

5.82 

1.803 

300.26 

0. 1089E 

03 

0.3064E 

02 

63 

6. 10 

i.  921 

300.2/ 

0.1136E 

03 

0. 3353E 

02 

64 

6.37 

z.OJi 

3  o  0 . 2  5 

0 . 1 2  8  OE 

03 

0. 3648E 

02 

65 

6.  70 

2.  173 

30U.25 

0. 139/E 

0  3 

0.4015E 

02 

66 

7.04 

2.326 

300. 24 

0. 1522E 

03 

0.4414E 

02 

67 

3.  U8 

1 . 0^6 

300.25 

0.32  73E 

02 

G.8381E 

01 

68 

3.64 

1  .  J.39 

300.26 

0 . 4  5  o  0  E 

02 

0.  U96E 

02 

69 

4  .  1  8 

1.253 

300.28 

0. 5803E 

02 

0. 1552E 

02 

7  C 

7.09 

2.34  7 

300.43 

0. 1321E 

03 

0.4408E 

02 

71 

3.07 

1.539 

3  uG . 47 

0. 8494E 

02 

0.2335E 

02 

72 

3.47 

1 .103 

300. 55 

0.4161E 

02 

0. 1082E 

02 

73 

2.96 

1 .035 

300.46 

0 .30226 

02 

0 • 7o64E 

01 

74 

2.68 

1 .005 

300.77 

0 .24616 

02 

0.6194E 

01 

75 

2 .48 

0.987 

300.65 

0.2089E 

02 

0. 5194E 

01 

76 

2.29 

0. 972 

300.61 

0. 1745E 

02 

0. 4280E 

01 

7  7 

7.53 

2.569 

300.27 

0 . 1 6  8  7£ 

03 

0.5006E 

02 

76 

8 .00 

2.  799 

300.36 

0 . 18  786 

03 

0.5623E 

02 

79 

8.51 

3.05  7 

300.41 

u. 20876 

03 

0.6309E 

02 

80 

2  1.  72 

12.019 

299.o2 

0 .89336 

03 

0.3272E 

03 

81 

22.2b 

12.495 

300.35 

0.92616 

03 

0. 3403E 

03 

8 

2  j  .  64 

13. 543 

300. 26 

0. 10066 

04 

0.3755E 

03 

83 

2  5.51 

15. 02  6 

299.84 

0. 1120E 

04 

0.4251E 

03 

84 

26.  b  1 

15. 979 

299. 42 

0.11966 

04 

0.4632E 

03 

o  5 

2  8.06 

16. 993 

300. 18 

0.12796 

04 

0.4988E 

03 

86 

2  9.60 

18.299 

299.90 

0.1  3  78E 

04 

0 . 5h50E 

0  3 

87 

30.17 

18. 733 

299 .34 

0. 14156 

04 

0.5587E 

03 

68 

14.84 

6. 894 

300.23 

0.50956 

03 

0.17276 

03 

.  H  -.1*. 
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0—16 

TAoLt  15 

EXPERIMENTAL  GAS 

,  FLOW  DATA  1 

FOR 

CORE  2 

6  l\  L »-  A  C  t 

OATA  FUk 

i> ANOfaTONE  SAMPL 

E 

kUI\ 

C6LP 

PMAk 

K  A 

D6LPM 

WAV 

(ATM) 

( ATM-1 ) 

(MD) ( ATM 

6EC/CM2  OP) (GM/SEO  CM2  CP) 

1 

0 • 74d5t 

00 

0 .  7  3566 

00 

0.4759E 

02 

0 . 2 1 0 1 E 

02 

0.91506-02 

2 

0.6322c 

u0 

0. 79256 

Oo 

0.47776 

02 

0. 20936 

02 

0. 7  392E-0  2 

3 

0.30916 

0  0 

0.63616 

00 

0.48696 

02 

0.2 045 E 

02 

0.37726-02 

4 

0 . 4224c 

60 

C. 8o936 

00 

0.4632E 

02 

0.2069E 

02 

0.45356-02 

5 

C . 395  1c 

00 

0.66106 

Ou 

0.49966 

02 

0. 2001E 

02 

0.4346E-02 

6 

0  . 3 5  3 4 1 

oO 

0 . 69946 

00 

0. 50036 

0  2 

0. 19996 

02 

0.38126-02 

7 

0.^2596 

uo 

0.91126 

Ou 

U. 49086 

02 

0.203  76 

02 

0.34026-02 

fa 

0.27956 

00 

0.9372E 

00 

0. 50676 

02 

0. 1973E 

02 

0.29336-02 

9 

0.23646 

00 

0.95976 

00 

0.50906 

02 

0. 19646 

02 

0.2391E-02 

10 

0 • 1 9o  06 

00 

0.96856 

CO 

0 . 5052E 

02 

0. 19796 

02 

0. 19446-02 

11 

0. 17186 

00 

0. 10016 

01 

U.5062L 

02 

9.19796 

02 

0. 16836-02 

12 

0 • o  0  4o  6 

ou 

0 . 7465E 

00 

0.47516 

02 

0.21056 

02 

0 • 99  3oE-02 

13 

0. 7623L 

00 

0 . 73956 

00 

0 .48226 

02 

0 . 2  0  74  E 

02 

0.94126-02 

14 

0 . 6  fa  4 1 1 

00 

0 .  7  o4b  c 

00 

0.47846 

02 

0 . 2090E 

02 

0.81096-02 

15 

C  •  3  o  fa  0  6 

0O 

U.8247E 

00 

0.4887E 

02 

0.2046E 

02 

0 • 65476-02 

lo 

0  •  4  6  5  0  E 

00 

0 . 66  366 

00 

0.49366 

02 

0. 2026E 

02 

0.51736-02 

17 

0 • 3 fa34t 

OU 

0. 69706 

GO 

0.49756 

02 

0.2010c 

02 

0.4139E-02 

lfa 

0  •  4  1  /oL 

uo 

u  .  6  fa  7  5  E 

ou 

0.4956E 

02 

0.2017c 

02 

0.45536-02 

19 

0.3364c 

00 

0  •  9  0  3  9  6 

00 

0.49506 

02 

U.2020E 

02 

0 . 3  784 E-02 

2U 

c  .  2  9  2  3 1 

uo 

0.93316 

00 

0.5037c 

02 

0.19786 

02 

0.3083E-02 

21 

0.32446 

01 

0*  33246 

UO 

0.40986 

02 

0.2440E 

02 

0. 7737E-01 

22 

0 • fa  1 1 5E 

00 

0  •  7  4  4  o  6 

GO 

0.47226 

02 

0.21 18E 

02 

0. 10046-01 

23 

0. 1212E 

01 

0.63866 

00 

0 • h  3  70E 

02 

0.2188E 

02 

0. 1690E-01 

24 

0. 1 39ot 

01 

0. 65906 

00 

0.44486 

02 

0. 2243E 

02 

0. 24756-01 

25 

0. 1  786  E 

01 

0.62546 

CO 

0.43906 

02 

0.22  78 E 

02 

0.29106-01 

26 

u . 14036 

01 
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TAbLE  15 

EXPERIMENTAL  GAS  FLOW  DATA  FOR  CORE  2 
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TA6LE  16 

EXPER  I MtNTAL  GAS  FLOW  DATA 

FOR  CORE  2 

LA  1  A  FUR 

SANDSTONE 

SAMPLE  FROM  SECTION  1 

RUN  P  1 A 

(  ATM  ) 

P2A 

(ATM) 

TMK  DO 

( K, )  (CC/SEC) 

DEL  P2 
( AT  M2 ) 

1 

10.93 

6.305 

300. 30 

0.3315E 

03 

0.  / 9  7  6  E 

02 

2 

9.28 

5.113 

299.84 

O.2506E 

03 

0. 5999F 

02 

3 

14.47 

6. 771 

299.92 

0.50  7 4E 

03 

0. L324E 

03 

A 

12.67 

7.410 

300.  16 

0.4116b 

03 

0. 1055E 

03 

5 

7  •  7  C 

4 .093 

299. 81 

0 . 1 8 i IE 

03 

0.4253E 

02 

6 

8.82 

4 .  <6  0  3 

299.40 

0.2161E 

03 

0. 8133 E 

02 

7 

14.  33 

8.835 

299.08 

0.510JL 

03 

0. 1334E 

03 

8 

lo  •  ^4 

10.115 

299.49 

0.6044b 

03 

0. 1614F 

03 

9 

1  /.89 

11.441 

299. 36 

0.6999b 

03 

0. 1891E 

03 

10 

19.16 

12.394 

299. 42 

0.7713E 

03 

0.2137b 

03 

11 

. .  20. 7  3  .. 

13.670 

. 299.51 

0.8656b 

03 

0.2436E 

03 

12 

2  3.2b 

1 3 . 7o2 

299. 22 

0. 1020E 

04 

0. 2936E 

03 

13 

24.61 

17.033 

299.74 

0.1115E 

04 

0. 3248E 

03 

14 

^o.  2  1 

18.  195 

299.46 

0. 1203E 

04 

0. 3558E 

03 

13 

2  7  *  64 

19. 420 

299. 12 

0.1296E 

04 

0. 3866E 

03 

16 

21.63 

19.330 

298. 88 

0.9183E 

03 

0.2627E 

03 
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r  abl  t 

16 

EXPER 1  MENTAL  GAS  PLOW  DATA  FOR  CORE  2 

uAi  A 

FUR 

S  Ai\uS  T  Ui\E  SAMPlE  FROM  SECTION  I 

RUN 

DcLP 

PM  AR 

KA 

OELPM  k‘AV 

(  ATM) 

i ATM-1 ) 

( MU) (ATM 

SEC/C  M2  C  P  )  ( GM/ SEC  CM2  CP) 

i 

0 . 4  6  2  / 1 

01 

0.1160E  00 

0 . 1 436  E 

02 

0 • 6  962  E 

02 

0.3559E 

00 

2 

0. 4io6t 

01 

0.1389E  U0 

0.1437 E 

02 

0.6957E 

02 

0.2696E 

00 

3 

0. 5699E 

01 

0*86  06  c~  0 1 

0. 1324E 

02 

0. 7 553 E 

02 

0.5437E 

00 

4 

0. 52  57E 

01 

0.9962E-01 

0. 1348E 

02 

0. 741SE 

02 

0.4414E 

00 

3 

C . jc 06E 

01 

0 • 1696c  00 

0. 1464E 

02 

0.6829E 

02 

0. 1952E 

00 

6 

0. 39 IjE 

01 

0 • 1 525E  00 

0. 1 44 5 E 

02 

0. 6922E 

02 

0.2329E 

00 

7 

C . 5659 E 

01 

0 • 8544E-  C  1 

0. 1314E 

02 

0. 7607E 

02 

0.54  76E 

00 

8 

0 . o 1 24E 

01 

0.  /369E-C1 

0.  1292E 

02 

0  .  7  7  3  6  E 

02 

U.6474E 

00 

9 

U  .  o44  7  E 

Cl 

0.O819E-G1 

0. 1279E 

02 

0. 78 i 7E 

02 

u .  7486E 

00 

10 

0  •  6  7  7  IE 

01 

0 . o3  3  7E  —  01 

0. 1248E 

02 

0.8015E 

02 

0.8243E 

00 

li 

0. 70 77t 

01 

0.5611E-01 

0. 1230E 

02 

0.6127E 

02 

0.923  7E 

00 

12 

0  •  7  5 1 9  E 

01 

0.5123E-01 

0. 1203E 

02 

0. 83 1 IE 

02 

0. 1087E 

01 

13 

0 •  7  737E 

01 

O.4777t-01 

U. 1 1 94E 

02 

0. 8374E 

02 

0. 1 185E 

01 

14 

G  .  c  0  1 2  E 

01 

0  •  h  d  04 1  —  01 

0.11 75E 

02 

0.8508E 

02 

0. 12786 

01 

15 

C.82  17 E 

01 

0.4250E-01 

0. 1 164E 

02 

0. 8593E 

02 

0. 1376E 

01 

16 

0.  7  2  9  6  E 

cl 

0. 3556E-U1 

0. 1206E 

02 

0. 6289E 

02 

0.9808E 

00 
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T  AbL  L 

1/  EXPERIMENTAL  GAS  FLLM  DATA  FOR  CORE  2 

DATA 

FUh  j ANOSTuNE  SAMPLE  FROM  SECTION  2 

RCN 

PI  A 

P2  A 

TMK 

QO 

DLLP2 

(  ATM  J 

(ATM) 

IK) 

(CC/SEC ) 

( ATM 2 ) 

L 

6*30 

6.402 

300.30 

0.3313E 

03 

0. 9480E 

,01 

2 

3.11 

4.392 

264. 84 

0.2306E 

03 

0. o83oE 

01 

3 

8.7  7 

7.82  8 

299.92 

0.5074E 

03 

0. 1364E 

02 

A 

7.41 

6 . 3  3  o 

300.  L  6 

0.4116E 

03 

0. i 193E 

02 

5 

4.09 

3.44  6 

299.81 

0 . 1  8  1 1 E 

03 

0.48  7  4 E 

01 

6 

4.60 

3.909 

299.40 

U.2161E 

03 

0. 5908F 

01 

7 

8.  do 

7 . 8  96 

299.08 

0.3100E 

03 

0. 1607E 

02 

8 

10.1  1 

3.183 

299.49 

0.6044E 

03 

0 . 1 88oE 

02 

9 

11.44 

10. 400 

299 . 36 

0.6999E 

0  3 

0.2274E 

02 

iO 

12.39 

1 1 . 326 

299.42 

0.7713E 

03 

0.2534E 

02 

11 

13.67 

12. 6  78 

299.31 

0 . 86  3qE 

03 

0.2867E 

02 

12 

13.  76 

14.383 

299.22 

0 . 1020E 

04 

0.33  72 E 

02 

13 

17.03 

13. 8l  7 

299. 74 

0.11156 

04 

0. 4071E 

02 

14 

18.19 

1  o  .  9  6  0 

299.46 

0. 1203E 

04 

0 • 4342  E 

02 

15 

19.42 

18. 183 

299.12 

0. 1296E 

04 

0 • 4644E 

02 

16 

14.33 

18.231 

298.88 

0.9183E 

03 

0. 3087E 

02 
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TAoLE 

1  7 

tXPtKiMtNTAL  GAS  FLJw  DATA  FOR  CORE  2 

^  1  M 

FOR 

8ANUSTLN6  5  AMPLE  FROM  SECTION  2 

ROlN. 

otLF 

PM  AR 

KA 

OfcLPM  WAV 

(ATM) 

(ATM-1 ) 

( MO) { ATM 

S6C/CM2  CP) (GM/SFC  C 

1 

0. 60296 

uO 

0.16946  00 

0.94776 

02 

0.  10336 

02 

0. 35o76 

00 

2 

0. 7213c 

00 

0.21046  00 

0. 96666 

02 

0.  1013E 

0? 

0.2702E 

00 

3 

0.94246 

00 

0.12036  00 

0.6767E 

02 

0. 1138E 

02 

0. 34346 

00 

4 

0.63406 

00 

0.14326  00 

0.9392c 

02 

0. 10696 

02 

0.4426E 

00 

3 

C • 0^-646 

00 

0 . 2o  3  3 6  GO 

0. 10036 

03 

0.99  726 

01 

0.19536 

00 

0 

0  .  o  9  4 1 6 

00 

0.23  306  GO 

0.98466 

02 

0. 10136 

02 

0.23346 

00 

7 

0.93946 

00 

0.11946  GO 

0.63946 

02 

0.11696 

02 

u. 34936 

00 

6 

o  .  47^9t 

00 

0 • 1 039 t  00 

G  •  6  6  6  6  c 

02 

0 . 1 1346 

02 

0 • 64956 

00 

5 

o  •  1 0  4  16 

01 

0.91376-G1 

0 • 6  33o  6 

02 

0. 1200E 

02 

0.75116 

00 

10 

o  •  i  0  6  6  6 

01 

0.6432c-Cl 

G. 62426 

02 

0.12136 

02 

0.82  736 

00 

11 

. . V.  10926 

JH . 

0.  /cd206-u1 

u .  6  1  6  9  E 

02 

0. 12216.. 

02  . 

0.92726 

OQ . 

12 

0.117/6 

01 

0.639 0c-0l 

0. 774JE 

02 

0.12916 

02 

0.10916 

01 

13 

0.1 2  36t 

01 

G  •  t>08  4  6  —  0 1 

0.74396 

02 

0.13416 

02 

0.11906 

01 

14 

0.12336 

01 

0 . 3o  69c- G i 

0.73396 

02 

0. i 32ot 

02 

Q.1283E 

01 

13 

L. 12336 

01 

0.33196-01 

0 . 7  36nt 

G2 

0.13196 

02 

0. 1382E 

01 

It. 

0.11196 

01 

0.72316-01 

o • 8  04 1 E 

02 

0. 1244E 

02 

0.96466 

00 
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TABLE 

IB 

EXPERIMENTAL  GAS  flow  DATA  FOR  CGRt  2 

DATA 

PuK 

SANDSTONE  SAMPLE  FROM  oeCTION  3 

K  U  N 

P  I A 

P  2  A 

TMK 

DO 

DELP2 

(  A  hvi ) 

(  ATM  ) 

(K) 

{ CC/bEC ) 

( ATM2 ) 

1 

3 .  :>0 

5.  14b 

300.30 

0.3315E 

03 

0. 3768E 

01 

2 

A  •  39 

4 . 093 

299.84 

0 • 2  506E 

03 

0. 2340E 

01 

3 

7.83 

7.476 

299.92 

0.3074E 

03 

0. 3364E 

01 

A 

6.36 

6.183 

300.  16 

0 . 4  L 16E 

03 

0.4725E 

01 

5 

3.43 

3.149 

249.81 

0. 1811E 

03 

0. 19o3E 

01 

6 

3.91 

3.399 

299. 40 

0 . 2  1 6 1 E 

03 

0.2323E 

Oi 

7 

7.87 

7.486 

299.08 

0.3100E 

03 

0 . 593  IE 

01 

8 

9.13 

8.720 

299.49 

0.6044E 

03 

0.741  IE 

01 

9 

10.40 

9.961 

299. 36 

0.6999E 

03 

0. 6937E 

01 

10 

11.33 

10. 880 

299.42 

0.7713E 

03 

0.9897E 

01 

11 

12. 38 

12.113 

299. 31 

0 . 86  3oE 

03 

0. 1 147E 

02 

12 

14.38 

14.  112 

299. 22 

0. 1020E 

04 

0. L35  7E 

02 

1-4 

13. 82 

13. 343 

299. 74 

0. 1 1 i3E 

04 

0. 14686 

02 

14 

16.56 

16.4/  / 

299 . Hb 

0. 1203E 

04 

0. 162  7E 

02 

13 

1  8.  1  8 

1  7  .  O  0  4 

299. 12 

0. 1296E 

04 

0. 1 794E 

02 

lo 

13.23 

12.777 

298.88 

0.9183E 

03 

0. 1181E 

02 
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EXPERIMENTAL  GAS  FLOW  OAT  A  FUR  CORE  2 

DATA 

FuK 

SANDS  TuNt  SAMPLE  FROM  SECT  IUN  3 

RUN 

DEEP 

PMAK 

KA 

DELPM  WAV 

(  ATM) 

(ATM-1 ) 

(MU) (ATM 

SEC/CM2  CP) (GM/SEC  CM2  CP) 

1 

U  •  3  3  3  6 1 

OU 

0.18786  00 

0.23836 

03 

0.4196E 

01 

0.3369E 

00 

2 

G  .  2  9  9  4  E 

00 

0.23376  CO 

u  •  2  6  o  2  E 

03 

0.  3  7  o6 F 

Q\ 

0.2703F 

00 

3 

0. 33C4E 

GO 

0.13G7E  00 

0 . 2  3o 1 6 

03 

0. 3904E 

01 

0.5457E 

00 

A 

u. 3709e 

UO 

Q.1370E  00 

0.2360E 

03 

0. 42386 

01 

0.4429E 

00 

_ _ & _ 

C  #.29  7  7 6 

00 

0.30336  00 

0.24886 

03 

0.40196 

01 

0.1936E 

00 

6 

0. 309oL 

00 

0  •  26646  00 

0.2302E 

03 

0.39976 

01 

0.2335E 

0  0 

7 

U  .  J  8  o  2  E 

00 

0 .  1 3 02 e  00 

0.2317E 

03 

0.43166 

01 

0. 5496E 

00 

6 

0 . 4 1 5 1 E 

00 

0.U2GE  00 

0.22046 

03 

0 • 4536E 

01 

0.6499E 

00 

9 

0  •  3  8  9  6 

00 

u.  9823E-01 

0.21206 

03 

0.4 71 76 

01 

0. 7317E 

00 

10 

0  *  4  4  5  7  fc 

00 

0. 90076-01 

0.2109E 

03 

0.47416 

01 

0. 82 79E 

00 

11 

U  .  4  t>  4  H  E 

00 

0. dl OOc- 01 

0 . 2096E 

03 

0.48876 

01 

0.92  786 

00 

12 

O.h  729e 

JO 

0 • 6  969  6—  0 i 

0.20376 

03 

0.49O9E 

01 

0. 1092E 

01 

13 

0-4712E 

UO 

0. 6418E-01 

u.2Go7E 

03 

0.4839E 

01 

0.U91E 

01 

14 

G . 4  8636 

00 

0. 396 IE- 01 

0  •  2  0  1  i  E 

03 

0.49746 

01 

0.12  84  6 

01 

IS 

C  *  3  0  (J  1 E 

00 

0.33  766-01 

0.  19626 

03 

0. 3097E 

01 

0. 13836 

01 

16 

0*43426 

GO 

0.  / o  90c— 01 

0.2 100E 

03 

0.47616 

01 

Q.9833E 

00 
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EXPERIMENTAL  GAS  FLG‘W  DATA  FOR  CORE  2 

u  AT  A 

Ilk 

SANDSTONE  bAMPLE  FROM  SECTION  4 

RUN 

Pi/. 

P2  A 

TMK 

go 

DELP2 

(  ATM  ) 

(  ATM  ) 

(KJ 

(CG/SEC ) 

( A TM2 ) 

1 

3.13 

4. 3b  1 

3o0. 30 

0.3315E 

03 

0 . 5b93E 

01 

2 

h.09 

3 . 506 

299.84 

0 . 2  5  OdE 

03 

0 . 4460 t 

01 

3 

7.48 

6  •  806 

299.92 

u. 3074E 

0  3 

0.9598E 

01 

4 

6.19 

3.381 

300.  16 

0.4116E 

03 

0. 7106E 

01 

5 

3.13 

L  •  664 

^99.81 

0 . 1  8  1 1 E 

03 

0.2818E 

01 

6 

3.6  0 

3.072 

299.40 

0 .2 16  IE 

03 

0. 35 18E 

01 

7 

/  .  49 

6*839 

299.08 

0.5100E 

03 

0 . 9260E 

01 

8 

6  •  {  2 

8.048 

299.49 

0 . 6  0  4  4  E 

03 

0. 1127E 

02 

9 

*.96 

9.272 

299. 56 

0.6999E 

03 

0. 1325E 

02 

10 

10.86 

10.191 

299.42 

0.771 3E 

0  3 

0. 1432E 

02 

11 

12.  11 

11.399 

299.31 

0 . 86doE 

0.5 

0. lbSOE 

02 

12 

14.  i  i 

13. 363 

299.22 

0. 1020E 

04 

0. 2037E 

02 

13 

13.33 

14.363 

299. 

0. 1 1 15E 

04 

0.2340E 

02 

14 

lb.  4d 

1 3. 6o6 

299.46 

0. 1203E 

04 

0.2544E 

02 

13 

17.68 

16.919 

299. 12 

0. 1296E 

04 

0.2O49E 

02 

16 

1^.78 

12.043 

298 . 68 

0.9183E 

03 

0. 1816E 

02 
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E  XP  Ek I MEN  T AL  GAS  FLOW  DATA  FOR  CORE  2 

DAI  A 

FOR 

SANDSTONE  SAMPLE  FROM  StCTiON  4 

RUN 

DcLP 

PMAR 

KA 

DtL  PM 

WAV 

(ATM) 

(ATM-1 ) 

( Mu) ( ATM 

SEC/CM2  CP ) ( GM/ SEC  CM2  CP) 

1 

U . 56o9b 

00 

0. 2  0  60 E 

00 

0. 1303b 

03 

0.  76  77E 

01 

0.3571c 

00 

2 

G • 5  Og 9E 

uU 

U.Z632L 

00 

0. 1253E 

03 

0. 7978E 

01 

0.2704E 

00 

3 

0.6720E 

Ou 

0* 1400E 

00 

U.1183E 

03 

0. 8453e 

01 

0.3439E 

00 

4 

G. 6039E 

00 

0. 17C0E 

00 

0. 1297E 

03 

0. 7712E 

01 

0.4431E 

00 

5 

6 . 4  a  4  8  E 

00 

0 • 34  4  1  E 

GO 

0 • 1*+ 33E 

03 

J.6979E 

01 

0.  19  d  7E 

a a 

6 

0.52  7 4E 

00 

0.2998E 

00 

0 . 1  3  6  6  E 

03 

0.7319E 

01 

0.23  3faE 

0  0 

7 

0 • 64o4E 

00 

0. 1396E 

GO 

o. 1226E 

03 

0 . 8  1  53  E 

01 

0. 5499E 

00 

8 

0.6720E 

uO 

0.  1193E 

GO 

0. 1 198E 

03 

J . 83h3E 

01 

0.6503E 

00 

9 

u • 6  fa  9  CE 

00 

0 • lu40c 

GO 

0.1  182E 

03 

0 . 84o3  E 

01 

0.7520E 

00 

10 

u . o69GE 

00 

u. 9492E- 

01 

0. 1 188E 

03 

0  .  fa  4 1 5  E 

01 

0.8283E 

00 

11 

u  .  7  i  4  3  E 

00 

0 . 83G6E- 

•01 

0 . 1 154E 

03 

0 . 8fao4E 

01 

G.9283E 

00  . 

12 

U.  7 485E 

Oo 

0.  72  79E- 

•Cl 

o.  ii  he 

03 

0.9003E 

01 

0. 1093E 

01 

13 

0 . 7  fa  2  5  fc 

00 

0 • Oo  8  7£- 

Cl 

0.10  71 1 

03 

0.9333E 

01 

0. 1 192E 

01 

14 

G.7910E 

00 

0 • 62 1 8  E- 

01 

0. 1062E 

03 

J. 9412E 

01 

0. 1285E 

01 

15 

C  .  /  6  5  3 1 

00 

0.5  780E- 

Cl 

0. 1G98E 

03 

0. 91 JbE 

01 

0. 1384E 

01 

16 

0. 7315E 

00 

0. 8057E- 

Cl 

0. 1 129E 

03 

0. 8856E 

01 

0.9859b 

00 
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1  Ad  LE  20 

EXPER I  MENTAL  GAS  F L 0 /w  DATA  FUR  CURE  3 

E  3  u  f  A  C  E  DATA 

FOk  L 1  ME  S I  ONE  SAMPLE  -  DRY 

RUN 

P 1 A 

P2  A 

TMK 

QU 

Ot  LP2 

(  ATM  i 

{  ATM  ) 

(K) 

(CC/SEC ) 

( AT  M2 ) 

i 

2  8.39 

6.  o92 

299.o4 

O.4670E 

0) 

0. 7609E 

03 

2 

61.29 

i.  /  •  3  3^. 

299. 37 

0 . 1  3  4  1 E 

04 

0.  i4  1  7 F 

04 

3 

3  7 .36 

16 . 332 

299. 05 

0. 1243E 

04 

0. 3041E 

04 

4 

34.07 

13. 307 

298.80 

0. 11 33E 

04 

0.2690E 

04 

3 

30.  74 

14,130 

298.93 

0. 1063E 

04 

0.2375E 

04 

6 

4  7.06 

13.046 

296.62 

0.9  796E 

03 

0. 2103E 

04 

7 

4  3.06 

11.673 

299.22 

U.873oE 

03 

0. 1770E 

04 

8 

3  9  .  o  3 

10. 187 

29d.83 

0. 7698 E 

03 

0.1 468F 

04 

9 

36*2  3 

9.214 

299.23 

0.6638E 

03 

0. 1229E 

04 

10 

32.33 

7.938 

299.40 

0 . 58596 

03 

0.9824E 

03 

11 

23.3  3 

3  •  1  9  o 

299.26 

0.3 744E 

03 

0. 5162E 

03 

12 

18.  3d 

3. 601 

299.47 

0.2681E 

03 

0. 3235E 

03 

U 

14.40 

2. 7  80 

299.56 

0  •  1 8  8  4E 

03 

0. 1997E 

03 

14 

1. 33 

0.930 

299.01 

0. 1433E 

01 

0. 9040F 

00 

IS 

3.81 

0.970 

299.o5 

0. 1724E 

02 

0. 1 36 1 E 

02 
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I  AdLt  20 

EXPERIMENTAL  OA5  FLOW  DATA  FOR  CORE  3 

t  ft  0 F  A  C  E 

DATA  FOR  LiMEbTUNE  SAMPLE  -  DRY 

RUN 

DEEP 

PMAR 

KA 

DELPM  WAV 

i  ATM) 

1 ATM-i ) 

l MO) ( ATM 

SEC/CM2  CP) (GM/SEC  CM2  CP) 

1 

0. 2169E 

02 

C. 5702E-01 

U.8103E 

01 

0. 1234E 

03 

0.5394E 

00 

2 

0 « 4342E 

02 

C.2327E-01 

0. 5  0 36 1 

01 

0.  1985E 

03 

0. 1434 E 

01 

3 

0.41 C2t 

U2 

o . 26  98L- u 1 

0.5272E 

01 

0. 1897E 

03 

0.  1358E 

01 

4 

0  .  3  b  7  /  £ 

02 

0. 2383E-01 

0.54  b6  t 

01 

0. 1823E 

03 

0.1259 E 

01 

5 

C .  36t>  It 

02 

0. 30b3E-Cl 

0.5  720E 

01 

0. 1748E 

03 

0.  Uoit 

01 

6 

0 . 3463E 

02 

0. 3294E-01 

0.5929E 

01 

J.  1687 1 

03 

0. 1074b 

01 

7 

0. 3 1 9  9  E 

02 

0.3614E-C1 

0 . 62  8  bt 

01 

0.  1 590  E 

03 

0.9593E 

00 

8 

0 • 2946t 

02 

0  .  'tO  1  3  E—  01 

0.6O49E 

01 

0. 1 304E 

03 

0.8484E 

00 

9 

0.2703L 

02 

0.4400E-01 

0.7059E 

ol 

0. 14 1 7  E 

03 

0. 7544E 

00 

10 

C.2439E 

02 

0.4966E-01 

0 . 7  5  5  7  E 

01 

0. 1323E 

03 

0.6477E 

00 

u 

0 . l b i 5E 

02 

0. 70 06 £—01 

0.9102E 

01 

0. 1099E 

03 

0.4163E 

00 

12 

0. 1458E 

02 

0.90 16E-01 

0.1 043  b 

02 

0. 9587E 

02 

0.2989E 

00 

13 

0.  llt>LE 

02 

0.1164E  00 

0. 1 185E 

02 

0. 8436E 

02 

0.2 1J4E 

00 

14 

L  •  4  J  U  uL 

00 

u  .  8  8  5  0  E  CO 

0. 1999E 

02 

0. 5002E 

02 

0. 1 597b- 

02 

15 

0 .2b45E 

01 

0.4180E  00 

0. 1 59b  E 

02 

0. 6262E 

02 

0. 1922E- 
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TA3L8  21 

8XP  tK I  MENTAL  bAS  FLOW  DATA 

FOR  CORE  3 

DATA  FOR  LIMESTONE  SAMPLE 

-  i-lOUlD 

SATURATION 

.0765 

RUN 

P 1 A 

P2  A 

TMK 

QO 

0ELP2 

(  ATM) 

l  ATM  ) 

(  K  ) 

(CC/S8C ) 

(ATM 2) 

1 

15.40 

2.781 

209. 28 

0.18678  03 

0.2321E  03 

2 

2  1.27 

4. 1 76 

2  99  9£1 

0.29448  08 

0.4349E  03 

3 

6  5  .  s  4 

i  7  .  74 

299.22 

0. 13298 

04 

0.39618 

04 

4 

62.  10 

16.499 

298.42 

0. 1252E 

04 

0. 33848 

04 

5 

5  3.71 

15.409 

298. 78 

0. 11638 

04 

0. 32098 

04 

6 

5  o  •  1 9 

14.273 

298.49 

0.10818 

04 

0.2842E 

04 

7 

51.82 

13.225 

298.66 

0.99808 

03 

0.25108 

04 

8 

46.6b 

11.591 

298. bl 

0.86818 

03 

0. 2042E 

04 

9 

4  2.46 

10.298 

298. 53 

0.76758 

03 

0.  16978 

04 

10 

3o  .  4  o 

9.074 

29o.95 

0 • 6  740 8 

03 

0.13998 

04 

11 

3  5.22 

8 . 10  7 

298.34 

U. 59938 

03 

0.11748 

04 

12 

31.20 

6.93  7 

299. 18 

0.50718 

03 

0. 9255E 

03 

13 

2  6  •  4H 

b.  383 

299.18 

0.4032E 

03 

0.66788 

03 

14 

1.48 

0.935 

299. 83 

0.17338 

01 

0. 1331E 

01 

lb 
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TAbLt 

21 

1XPER I  MENTAL  GAS  FLOW  DATA  FUh  CURE  3 

DATA  F(JK 

LlMESTONt 

SAMPLE  -  LiwllIO  SATUKAT  IUN  .0765 

HUM 

DLL*- 

PM  AH 

KA 

DEL PM  WAV 

(  AT  h ) 

(ATM-1 ) 

<  MD) ( ATM 

SEC/CM2  CP) (GM/SEC  CM2  CP) 

1 

0. 12  7Jt 

02 

0.10951  00 

0. 10101 

02 

0.99021 

02 

0.2086E 

00 

2 

C.l 7C91 

02 

0.  76601-01 

0.85321 

01 

0.  1 1 72h 

03 

0.32  77F 

00 

3 

0.47 77E 

02 

0 .2412E-C1 

0.43351 

01 

0.23071 

03 

0. 1439E 

01 

4 

0.43601 

u2 

0.2545E-01 

u .  448 1 1 

01 

0.22311 

03 

0.  1 36  IE 

01 

JL 

0.43  3.0  fc. 

QZ 

0.2 6 9oi-01 

0 . 4o6bL 

01 

0.21421. 

.03 

0.  1 2  721 

01 

6 

0.40921 

02 

O.28791-01 

0.48351 

01 

0.20591 

03 

0.1180E 

01 

7 

0. 36601 

02 

0. 3075c-01 

0 . 5  0  7 1 1 

01 

0. 1972E 

03 

0. 1092E 

01 

b 

C . 35061 

v2 

0. 3434E-01 

0. 54011 

01 

0. 1852E 

03 

0.9535E 

00 

9 

0.32171 

62 

0.37911-61 

0.57311 

01 

0. 1745E 

03 

0.8453E 

00 

10 

0.29411 

02 

0 . 42  Cot- 01 

O.olObl 

01 

0.  16371 

03 

0. 74336 

00 

LI 

_  ,.0.2/111 

02 

0 . 4  u lot-01 

0 . 64321 

. 01 

0. 15501 

03 

0.66251 

00 

i2 

0.24271 

02 

0.32  44 E—  01 

0.69281 

01 

0. 1443E 

03 

0.3615E 

00 

13 

0 . 2  Gout 

u2 

0 • o24ot— 01 

0. 76131 

01 

0. 13141 

03 

0.44771 

00 

14 

0.53001 

00 

0.62641  CO 

0 . 16411 

02 

0. 6090E 

0? 

0. 19321- 

•02 

15 

C • 40991 

01 

0 • 32  4oE  CO 

0. 132bE 

02 
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TABLE  22 

EXPERIMENTAL  bAS  FLOW  DATA  FUR  CURE  3 

LAI m  FUR  LIMESTONE 

SAMPLE  -  LlOUiD  SATURATION  ,124 

RUN 

P  1 A 

P2  A 

TMk 

QO 

DELP2 

l  ATM) 

l  ATM  ) 

IK) 

tCC/SEC) 

( A  T  M2  ) 

i 

1  / .  3b 

2.838 

299.66 

0. 1938E 

03 

0. 2929E 

03 

2 

2  2.92 

3.994 

299.99 

0. 2d20E 

0  3 

0 . 4  8  6  7  F 

0  3 

3 

71.  96 

17.644 

300. 39 

0. 1338E 

0^ 

0. 4867E 

04 

4 

o 

• 

O 

1 6 • ^36 

299.93 

0 . 128  IE 

04 

0.4483E 

04 

5 

o  3  •  2  6 

13. o48 

2  9  9 . 6  7 

0. 1196P 

04 

0.4007E 

04 

6 

o2  *  1  o 

14.963 

299.26 

0.112/E 

04 

0.3640E 

04 

7 

3d. 69 

13.969 

296.94 

0. 1048E 

04 

0. 3249E 

04 

a 

34.  34 

12.792 

299.23 

0.9340E 

03 

0. 28 1  IE 

04 

9 

30.12 

11.333 

299.21 

0 • ob68E 

03 

0.2379E 

04 

ID 

‘♦3.13 

10. 138 

299.03 

0.7491c 

03 

0. 1936E 

04 

JU 

.  40.32  . 

6.771 

. 299.09 . 

0  •  6446L 

03 

0. 1 34  9 1 

04 

12 

3  3 . 0  ci 

7.342 

299.63 

0.3361c 

03 

0. 1177E 

04 

13 

31.47 

6.3  79 

299.63 

0.4626E 

03 

0. 9497E 

03 

14 

26.3  7 

5 . 04  6 

299.86 

0. 3616F 

03 

0.6697F 

03 

IS 

19-78 

3.470 

299.90 

0 .2415E 

03 

0.3794E 

03 

16 

1.  /  4 

0.936 

300. 23 

0.23  70t 

01 

0.2162c 

01 

17 

3.  94 

1.039 

300. 30 

0 . 28 70c 

02 

0. 3423E 

02 
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TAbL t  22  EXPERIMENTAL  GAS  FLOW  DATA  FOR  CURE  3 


DATA  Fuk  LIMoSTUNE  SAMPLE  -  LIQUID  SATURATION  .124 


RUN  DELP  PM AR  KA  OELPM  WAV 

(ATM)  (ATM-1)  ( M D )  (ATM  SEC/ CM2  CP)(GM/SEC  lM 2  CP) 


1 

U  .  i.4 49E 

02 

0.98960-01 

0.8  339E 

01 

0. 1199E 

03 

0.2160E 

00 

z 

0 . 1  b4;>E 

02 

0. 7d71E-01 

0 . 7  3  2  4 1 

01 

0. 1 3o5E 

03 

C. 31 33E 

00 

3 

0. 3431E 

U2 

0 . 2232  E~  0 1 

0.35900 

01 

0.27860 

03 

0. 1439E 

01 

4 

C. 5213b 

02 

u.232oE-0i 

0.3704E 

01 

0.2700E 

03 

0. 1384E 

01 

C .  4  9  4  1 L 

02 

0. 24O60  —  01 

0.38640 

01 

0.2588E 

03 

0. 1295E 

01 

6 

0.4  720E 

02 

0. 25930-01 

0.3988E 

01 

0. 2507  E 

03 

0.1224E 

01 

7 

0.44  72E 

02 

0.27530-01 

0.4 1370 

01 

0.2417E 

03 

0.11 42 E 

0  1 

8 

0.41  /5t 

02 

0.2970E-01 

0.43490 

01 

0.2299E 

63 

0. 10410 

01 

9 

0. 3638E 

02 

0. 3244E-01 

C • 4  6  03  E 

01 

0. 2172o 

03 

0.93  79 E 

00 

10 

0  •  o  b  0 1 E 

02 

0 . 3o  1  8E-C1 

0.4925E 

01 

0.2030E 

03 

0.8229E 

00 

11 

o  .  3  1  630 

U2 

0  •  4  0  74  E—  0 1 

0.32850 

01 

0. 1892E 

03 

0. 71020 

00 

12 

0.2774 E 

02 

0.47130-01 

0 . 5  7  8  7  0 

01 

0.17280 

03 

0.3917E 

00 

13 

G.23C9E 

02 

G.5284E-01 

0.61810 

01 

0. 1 6 1 8E 

03 

0.5115E 

00 

14 

U. ,01320 

02 

O.u36bh-Oi 

0 . 66380 

Oi 

0. 1462E 

03 

0.4012E 

00 

15 

0  . 

02 

0.8O01E-01 

0.8035E 

01 

0. 1245E 

03 

0.2688E 

00 

16 

0 . 3  U  7  0  E 

00 

0 •  7465 1  00 

0. 1382E 

02 

0. 7233E 

02 

0.26430- 

■02 

1  7 

0.4  9  03E 

01 

0.2865E  00 

0. 1 0590 

02 

0. 9455  E 

02 
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fcXPERIMENTAL  GAS  FLUrt  U  AT  A  FOR  COKE  3 

LATA  FuR 

L i KESTGNL 

SAMPLE  -  LIUUID  SATURATION  .203 

RUN 

Pi  A 
(  ATM  ) 

P2  A 
{ATM  ) 

TMn 

*3 

ICC/SEC) 

0ELP2 

( A  TM2 ) 

1 

79.43 

17.  174 

299. 3b 

0 . 1 3  64t 

04 

0.6017E 

04 

2 

7  9.30 

10.209 

298.37 

0.1291E 

04 

0.3436E 

04 

3 

7^.30 

13.486 

295.61 

0 . 1227t 

04 

0.4988E 

04 

A 

6  7.  d4 

14.261 

235. 70 

0. 1 L28E 

04 

0.43586 

04 

3 

61.83 

12.519 

298,39 

0. 1009E 

04 

0. 3o58E 

04 

6 

3  7. 03 

11.387 

298.54 

0 . 90o  7 1 

03 

0. 31 18E 

04 

7 

3  1 . 9o 

10.313 

295. 56 

0. 8049E 

03 

0.2393E 

04 

8 

4 1  •  66 

9.042 

298. 86 

0.7000E 

03 

0.2114E 

04 

9 

4  1  •  o2 

7.  763 

2  99.  19 

0 . 39  77E 

03 

0. 1672E 

04 

10 

3  o  •  b  2 

0.627 

299. 34 

0. 3039E 

03 

0.1312E 

04 

ii 

5  1  •  5  b 

3.443 

299. 24 

0.4106E 

03 

0.9741E 

03 

12 

2o •  3  7 

4.279 

299. 36 

0.31 76L 

03 

0.6773E 

03 

13 

19.64 

A.  932 

299. 32 

0.2077E 

03 

0.3770E 

03 

14 

Z  .  jO 

0.940 

299.81 

0.2776E 

01 

0. 3108E 

01 

13 

6.93 

1.062 

299.53 

0 • 3066E 

02 

0.4177E 

02 

i  J  jA  T 


J  -  >‘1AU  .a^  J  ^  (d  AO  JAT.i;;  I,-.JSAj 


c  A. 

- 

a  In  J 

i< 

(  r.  ) 

A  IS 

t  >  3  c  \  J  i  l  /U  l  I  m  )  (  |\j 

AU 

•  \  i  /.  . 

t4  0 

.1  i  . 

;  ^  .  c 

A  \  i  .  \  i 

d  t-  .  \ 

i 

A'. 

-  . 

fcG 

.  ;  ■  i  . 

\  C  .  ji'i 

a  t  .  J 

'MU 

£  , 

•  C 

+  . 

3  '■  i  .  U 

lo.  ^ 

f  »  .  v!  i 

Qfc 

*  0 

t  C1  1  4  . 

AO 

.<  11. 

)  \  •  o  d  S 

1  i  :»  •  •*  1 

AC  .  \  0 

A 

ct  i  . 

J  . 

'  .  V 

i  •  .a 

vO  •  1 0 

AC 

til 

J  1  ■ . .  •  •  u 

A  •  •  .  \ 

\  .  i 

t  U  .  \  d 

A 

1  c  .  .  C 

J ' 

OC  .  t  .> 

t  1  .  u  i 

ur'  .  1  c 

\ 

■M 

-  I  i  .  . 

\  ,0 

.  ■ 

uo.  a 

■>  O  .  0  A 

Al 

N  .  i  . 

^  \  <  . 

. 

d  a  \  .  v 

i>  a  •  i  a 

>  0 

i  •  1  <  a  . 

'  c.  <  .  0 

At  .  v  a 

t  ^ •  o 

• 

01 

1m  .  „ 

c  L-  1  •  .  W 

• 

t  i  .  c 

.  .  e 

1  i 

i  <  '•  ^  y  . 

\  1  t  .  V< 

A  \  .  t 

\  c  .  oS, 

jCf  U  .0 

\  \  A  .  U 

s>  . 

.  <.  .  ^ 

A  U  .  t  1 

ti 

It 

(.  i  <  .  J 

1 

. 

1  .  v  . 

ww  .  ^ 

Ai 

•  \  ^  ;  ♦  . 

.•  «  .  a 

:  .  '  _> 

•  J 

cc  #d 

e  i 

L)—  3  3 


T  AbLh 

23 

EXPERIMENTAL  GAS  FLUw’  DATA  FOR  CORE  3 

DAT  A  6uK 

LIMES!  ONE 

SAMPLE  -  LIQUID  SATURATION  .203 

RUN 

DtLP 

PMAR 

KA 

DELPM 

«AV 

(  ATM) 

(  ATM-1 ) 

i MO) (ATM 

S6C/CM2  * 

CP) (GM/SEC  CM2  CP) 

1 

0.O228E 

02 

O.2070E-01 

J.2954E 

01 

0.33856 

03 

0. 1 4o4E 

01 

2 

0.;92JE 

02 

0.21  79E-01 

0. jG  /  it 

01 

0. 3256E 

03 

0.13936 

01 

3 

u . 3o82t 

02 

0.2278E-01 

0.31 75E 

01 

0. 3130E 

03 

0.13276 

01 

4 

0  •  3  3  2  8  £ 

02 

0. 2443c— 01 

0.3330L 

01 

0 . 30u3  c 

03 

0.12236 

01 

3 

L  .  4  0  0  i  c 

02 

0  #  2  6  7  i£~  0 1 

0.3330F 

. 0.1 

Q .2  8  3  3  £ 

03 

0. 1099E 

01 

t> 

0  .  *+  3  4  4  L 

02 

0.2915E-01 

0.3728 E 

01 

0.2682E 

03 

0 • 992  1 E 

00 

7 

0.410^0 

02 

0.32 126-01 

U.3952E 

Cl 

0. 25306 

03 

0 . 8  8  2 1  £ 

00 

6 

0.3781E 

02 

u.  33  78E-01 

0.42126 

01 

0.23746 

03 

0. 76916 

00 

9 

0.3385L 

02 

0. 4030 E- 01 

0.45436 

01 

0.2201E 

03 

0.6  5836 

00 

10 

0 . 30196 

02 

0 . 46U3L-01 

U  .  4 o  92  6 

OL 

0.20446 

03 

0.5585E 

00 

JJL 

0. .  2  6.2  4  c 

02 . 

0.336 7 t-Oi 

0.3328L 

.01 

0 .18  77 1 

03 

.0.43.4  8  E 

00 

12 

0.22096 

02 

0. 6  3  2  3  E  —  01 

0.5917E 

01 

0. 16906 

03 

0.3527E 

00 

13 

0 . 1 6  7 1 E 

02 

0. 8862E-01 

U. 69286 

01 

0. 1443E 

03 

0.23  166 

00 

IA 

o . i U  58c 

ILL 

0.3d 07c  .00 

.  0*  112166 

0.2 

0 .  68  76E 

02 

0.3Q96E- 

lQ2 

0.3't886  Cl  0.2627c  CO  G.9250E  01  0.10816  03  0.3418E-01 
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EXPtRI MtNTAL  GAS  FLOW  DATA  FOR  CORE  3 

t,  A  f  A  PGR 

Lif'coTGNb 

SAMPLE  -  LIQUID  SATURATION  .233 

RUN 

Pi  A 

P2  A 

TMK 

go 

DELP2 

{  ATM) 

(AIM) 

{  K » 

(CC/SEC) 

( ATM2 ) 

I 

8  5.44 

15. 22  o 

2  68 . 29 

0 . 1209E 

04 

0. 7067b 

04 

2 

b  1  •  4  6 

14.444 

237.89 

0.1146b 

04 

0.6426E 

04 

3 

7o .  0 5 

13. 328 

298.28 

0 . 10  52E 

04 

0. 5o05E 

04 

4 

70.33 

12.218 

298.19 

0.9393E 

03 

0.4826E 

04 

3 

6  6.57 

ll.lo4 

297.90 

0.6  / 30E 

03 

0.4174E 

04 

6 

6i.U4 

10.259 

298.32 

0. 79  76E 

03 

0.3621E 

04 

7 

6  6  •  6  0 

9.163 

29b . 5b 

0. 7068E 

03 

0.3007E 

04 

8 

30.  19 

8.081 

298.24 

0.6223b 

03 

0.2454E 

04 

9 

‘tS.  42 

7.115 

298. 42 

0.5428b 

03 

0.2013E 

04 

10 

39.98 

6  •  G*+  7 

299.06 

0 . 4564b 

03 

0.  1 562  E 

04 

u 

3i  .08 . 

. 5. 1 15. 

. 2  99.03 

0 .3813b 

03 

0. 1203E 

04 

12 

29.81 

4.  141 

298.85 

0 . 30  36t 

03 

0 . 8  7 1 aE 

03 

13 

2.ob 

0.943 

299.06 

0. 3930b 

01 

0.6149E 

01 

14 

dj.u  J 
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299.01 
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TaDlE  24  EXPER  I  MENTAL  G  AS  FLOW  DATA  FOR  CORE  3 


DaIA  FOR  LIMESTONE  SAMPLE  -  LIQUID  SATURATION  .253 


RUN 

DEL  P 

PMAR 

KA 

DEL PM  WAV 

i  ATM) 

( ATM-1 ) 

( MD) i AIM 

SEC/CM2  CP ) ( GM/ SEC  CM2  CP) 

I 

0  .  I  0  2 1 E 

U2 

0 . 1 9  8  7  E  -  0  i 

0.2220E 

01 

0.4505E 

03 

0. 1296E 

01 

2 

0  .  o  7  0  i  c 

U2 

0.2  0  86E—  C 1 

0 .2303E 

01 

J.4342E 

03 

0. 1236E 

01 

3 

0  .  o  2  7  2 1 

02 

0.2238E-01 

0.2420E 

01 

0  •  4 1  8  3  E 

03 

0.1133E 

01 

4 

0.5832E 

u2 

0. 2417 t- Cl 

0 . 2  383E 

01 

0.3917b 

03 

0. 1042E 

01 

5 

0. 343ab 

02 

0. 2606E-Q1 

0 . 2o  77E 

01 

0.3  735E 

03 

0.9534L 

00 

6 

u. 3078E 

o2 

0. 2805E-G1 

0.2813E 

01 

0.3555E 

03 

0.8707E 

00 

7 

0 • 4o43L 

02 

0.3088t-01 

0.2996E 

01 

0. 3337 E 

03 

0.7737E 

00 

8 

0  »  4  2  1 1  E 

u2 

0.3432E-01 

0.3217E 

01 

0. 3109E 

03 

0  •  6  8  4 0  E 

00 

9 

o.  J831<_ 

02 

0. 33C7E-01 

0.3415b 

01 

0 . 292  BE 

03 

0.598  IE 

00 

10 

C.  3  J 93E 

U2 

U.4348E-01 

0.3  704 E 

01 

0. 2700E 

03 

0.5036E 

00 

11 

0.299 7  E 

02 

0.4976L-01 

0 .400 1 E 

01 

0.2499E 

03 

0.42  19c 

00 

12 

U.2567E 

02 

0. 5891 E- 01 

0 . 4  3  8  8  E 

01 

0.22  79 E 

03 

0.33  71E 

00 

13 

0 . 1  7  1 0  c 

01 

0.3362L  00 

0.8100E 

01 

0. 1234E 

03 

0.4404E- 

■02 

14 

0  •  7  3  32c. 

01 

U . 2067c  00 

o • u  96  7  L 

01 

0. 1433E 

03 

0.4495E- 
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